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We review a number of topics related to block variable renormalisation 
^ ' group transformations of quantum fields on the lattice, and to the emer- 

c3 ■ ging perfect lattice actions. We first illustrate this procedure by consid- 

ering scalar fields. Then we proceed to lattice fermions, where we discuss 
perfect actions for free fields, for the Gross-Neveu model and for a super- 
symmetric spin model. We also consider the extension to perfect lattice 
perturbation theory, in particular regarding the axial anomaly and the 
quark gluon vertex function. Next we deal with properties and applica- 
tions of truncated perfect fermions, and their chiral correction by means 
of the overlap formula. This yields a formulation of lattice fermions, 
which combines exact chiral symmetry with an optimisation of further 
essential properties. We summarise simulation results for these so-called 
overlap-hypercube fermions in the two-flavour Schwinger model and in 
quenched QCD. In the latter framework we establish a link to Chiral 
Perturbation Theory, both, in the p-regime and in the e-regime. In par- 
ticular we present an evaluation of the leading Low Energy Constants 
of the chiral Lagrangian — the chiral condensate and the pion decay 
constant — from QCD simulations with extremely light quarks. 
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Motivation and Overview 



Over the recent decades quantum field theory has been estabhshed 
as the appropriate formalism for particle physics, as far as it is explored 
experimentally. Its treatment by perturbation theory led to successful 
results, for instance in Quantum Electrodynamics (QED), in the elec- 
troweak sector of the Standard Model and in Quantum Chromodynamics 
(QCD) at high energy. However, there are still many open questions, 
which require results at finite coupling strength — beyond the range of 
perturbation theory — such as numerous aspects of QCD at low and 
moderate energy. 

A method is known which has the potential to provide fully non- 
perturbative results for a number of field theoretic questions. This me- 
thod applies Monte Carlo simulations to lattice regularised quantum field 
theories. The generic uncertainty of perturbation theory — uncontrolled 
contributions beyond the calculated order — disappears in this approach. 
However, one has to deal with statistical errors, as well as ambiguities in 
the extrapolation to the continuum and to a large volume. 

Simulation resTiIts are obtained at finite lattice spacing, which causes 
systematic artifacts in the numerically measured observables. The stabil- 
ity of dimensionless ratios of observables under the variation of the lattice 
spacing is denoted as the scaling behaviour. Its quality, which is vital 
for the reliability of the continuum extrapolation, depends on the way 
in which the lattice regularisation is implemented. This work deals with 
renormalisation group techniques to improve the scaling behaviour com- 
pared to the standard lattice formulations, which describes, for instance, 
derivatives simply by differences between nearest neighbour lattice sites. 
In contrast to Symanzik's program, this technique does not attempt to 
correct a specific order in the lattice spacing, but it directly addresses a 
finite cutoff. We describe the renormalisation group approach in detail 
and present a variety of results that it led to, in particular for fermionic 
systems. 

The symmetries of a model under consideration are a key aspect for 
a controlled continuum extrapolation. A great virtue of the lattice regu- 
larisation is the conservation of exact gauge symmetries. But global sym- 
metries are often explicitly broken by the lattice structure, for instance 
the continuous Poincare invariance. The question, how well — and if 
— they are restored as we approach the continuum limit is a notorious 
issue, which is related to the scaling behaviour. Again it depends on 
the features of the lattice formulation, i.e. on the extent of the explicit 
symmetry breaking due to a finite lattice spacing. The renormalisation 
group technique provides a tool to improve the symmetries on the regu- 
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larised level — in principle they can even be implemented exactly — 
which renders the continuum limit smoother and safer. This property is 
particularly relevant for the (approximate) chiral symmetry of (almost) 
massless fermions. The chiral symmetry is essential for instance in QCD 
at low energy, and its discussion will take a central role in this work. 
Here we also need a chiral extrapolation, in addition to the limits that 
we mentioned already. 

As an introduction, we summarise in Section 1 some basic aspects of 
quantum field theory in the functional integral formulation. In particular 
we sketch the road from classical mechanics to Euclidean quantum field 
theory, with quantum mechanics and classical field theory as intermediate 
steps. We focus on the lattice regularisation, which we first introduce for 
the case of scalar fields. This allows us to summarise the notions and 
notations used throughout this work. 

In Section 2 we describe — still for scalar fields — the concept of 
block variable renormalisation group transformations. Under iteration 
they lead to a perfect lattice action, which is free of any cutoff artifacts. 
We also encounter approximations to a perfect action, which are needed 
for practical purposes, such as the classical perfection and the truncation 
of the couphngs. Still the lattice artifacts can be kept small for such ap- 
proximations, as we illustrate for the dispersion relation, the topological 
susceptibility of a quantum rotor and for thermodynamic quantities. 

In Section 3 we proceed to fermionic quantum field theories, where we 
start again with a few generalities. We review the traditional formula- 
tions of lattice fermions and describe the doubhng problem. It is related 
to the difficulty to keep track of the chiral symmetry in a regularised sys- 
tem — an obstacle, which obstructs other regularisation schemes as well. 
Therefore the existence of light quarks is an amazing feature of Nature, 
which cannot be described easily in a natural way. In that context, we 
discuss a brane world scenario as a possible solution to this hierarchy 
problem. 

Section 4 applies the concept of perfect actions to lattice fermions, 
which takes us to the main topic of this work. For free fermions, we 
demonstrate that this approach provides both, a solution to the doubhng 
problem and at the same time an exact scaling behaviour. Depending 
on the choice of the renormalisation group transformation, we can im- 
plement either locality or standard chirality in the perfect lattice action. 
However, if we insist on locality, the resulting Dirac operator still obeys 
the Ginsparg- Wilson relation, which guarantees a lattice modified but 
exact chiral symmetry. 

In principle, exact scaling and chirality can also be realised at finite 



3 



lattice spacing in the interacting case, where, however, perfect actions can 
in general not be constructed explicitly. An exception is the Gross-Neveu 
model in the limit of a large number of flavours. Here we present a perfect 
action for staggered fermions, and we approve the perfect scaling for the 
ratio of the chiral condensate and the dynamically generated fermion 
mass. The concept we are dealing with also reaches out to perfect lattice 
currents. With that ingredient, perfect actions can even capture exact 
supersymmetry on the lattice. 

In Section 5 we consider perfect lattice perturbation theory. We give 
results for the anharmonic oscillator and the Yukawa term, which couples 
fermions to a scalar field. In lattice gauge theory we show that the perfect 
lattice action retrieves correctly the axial anomaly, and we discuss the 
quark gluon vertex function in QCD. 

For practical applications, i.e. for the applicability in simulations, the 
couplings have to be truncated. In Section 6 we describe our trunca- 
tion scheme for the perfect fermion to a so-called hypercube fermion, 
which has been simulated successfully in the Schwinger model. In QCD 
it has been used for the spectral functions at flnite temperature, and — 
together with a truncated perfect vertex function — in the evaluation 
of the charmonium spectrum. For truncated perfect fermions, the scal- 
ing behaviour and chirality are not exact anymore, but the latter can 
be corrected again by inserting the hypercube fermion into the overlap 
formula. 

This procedure yields the "overlap hypercube fermion", which is an 
exact solution to the Ginsparg- Wilson relation. Its construction and 
properties are presented in Section 7. Similarly we can arrange for a 
modified but exact parity symmetry for lattice fermions in three dimen- 
sions. In two dimensions we review simulations results for overlap hy- 
percube fermions in the two-flavour Schwinger model, which reveal an 
excellent scaling behaviour. Here and also in QCD we further observe a 
strongly improved level of locality and approximate rotation symmetry 
compared to the standard overlap fermion. 

Section 8 flnally presents simulation results with Ginsparg- Wilson 
fermions in QCD, using the overlap hypercube fermion as well as the 
standard formulation of overlap fermions, both in the quenched approx- 
imation. This enables simulations near the chiral limit. Here our main 
goal is a connection to Chiral Perturbation Theory. This is an effective 
theory of strong interactions at low energy, which provides a variety of 
successful predictions. However, its effective Lagrangian involves free pa- 
rameters denoted as the Low Energy Constants, which play an important 
role in the physics of light mesons. Their theoretical determination can 
only emerge from QCD as the fundamental theory. This is a challenge 
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for lattice simulations, and the principal issue of Section 8. 

Wc measured light meson masses in the p-regime (characterised by a 
large volume), and we reveal the difficulties to evaluate Low Energy Con- 
stants in that setting. Then we focus our interest on the e-regime, which 
deals with a small volume. In the e-regime, the topological sectors play 
an extraordinary role. Hence we first give results for the distribution of 
topological charges and the resulting susceptibility, which is relevant for 
the mass of the Vj meson. Next we describe a 3-loop calculation which 
confirms the perturbative renormalisability of the e-expanded effective 
theory. We then apply various techniques to extract the leading Low En- 
ergy Constants: the chiral condensate — which is the order parameter 
of chiral symmetry breaking — and the phenomenologically known pion 
decay constant. In particular, the density of low lying eigenvalues of the 
Dirac operator is fitted to predictions by chiral Random Matrix Theory. 
The axial-vector current correlator, as well as the zero-mode contribu- 
tions to the pseudoscalar density correlation, are confronted with for- 
mulae of quenched Chiral Perturbation Theory. We will see that these 
methods do have the potential to evaluate the Low Energy Constants 
with their phenomenological values — which correspond to the large vol- 
ume limit — even in the e-regime. However, the final results have to 
await the feasibility of dynamical QCD simulations with chiral quarks. 

Section 9 is dedicated to concluding remarks, summarising the status 
of the fields of research that we addressed, along with an outlook on 
future perspectives. 
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1 Introduction 



1.1 From classical mechanics to quantum mechanics 

In classical mechanics, the trajectory x{t) of a point particle between 
fixed endpoints x(0) and x{T) is — in simple situations — determined 
by the principle of least action, which imposes the condition 6S = 0. The 
action S[x] is a functional of the conceivable particle paths x{t), 

S[x] = [ dtL{x,S) , (1.1) 
where L is the Lagrange function. A simple form of it reads 

L{x,x) = -x{t)'-Vm) , (1-2) 

with the particle mass m and a potential V (which we assume to be 
velocity independent). The variational condition SS = corresponds to 
Newton's equation of motion, mx = — VV^, at each instant t G [0,T]. 

Let us consider this transition in quantum mechanics. In contrast to 
classical mechanics, we now deal with a transition amplitude, which picks 
up contributions from all possible paths connecting the fixed endpoints. 
Hence the path in between is not determined. These contributions are 
summed up coherently, 

(f(T)|x(0)) = J'^x exp [^S[x]^ . (1.3) 

This expression represents a path integral (or functional integral), where 
the functional measure T>x symbolises the summation over all possible 
paths (which formally requires an infinite dimensional integral) [1]. 

In the (hypothetical) limit /i — > solely the classical path (which we 
assume to be unique) contributes, whereas the additional contributions 
for h > correspond to the quantum effects. However, if a path far from 
the classical one is varied, the phase in eq. (11.31) tends to rotate rapidly, 
so that such contributions almost cancel. As long as h is small compared 
to the action shift caused by path variations on the scale of interest, it is 
the vicinity of the classical path that dominates the transition amplitude 

m- 

This situation has a historically older counterpart in optics, where 
the classical and the quantum mechanical description correspond to the 
principles by Fermat and by Huygens, respectively. 

In order to attribute an explicit meaning to the functional measure 
Vx, we divide the period T into N equidistant intervals of length a = 
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T /N . In this discretised system, the path integral is given by — 1 
integrals over the possible positions at the times tj = j -a, j = 1 . . . N — 1. 
The expression f ll.Sp is then understood as the continuum extrapolation 
a —>■ (which, at fixed T, corresponds to — *• cxo). 



N-l 

n 



d'xj ... ".Z I Vx ... (1.4) 

IR3 



1.2 Classical field theory 

In field theory we do not consider particle paths x(t), but instead fields 
0(a;), where x = (t,x) is a point in space-time. The (classical) field takes 
its value in some abstract space, like IR" or C", for example. Now space 
and time are treated on an equal footing (up to the signature in the 
metrics) , which is a prerequisite for covariance. Moreover, the number of 
degrees of freedom is extended drastically: before there were just three 
of them (in each time point t), but now there is a degree of freedom for 
each field component in each single space-time point x. 

We assume in each point x a Lagrange density £(0, d^cp) to be defined 
(yU = 0, . . . , 3), which we denote as the Lagrangian. The field theoretic 
action is given by 

^[0] = J d^xC{(P,d^(P) . (1.5) 

Now an action value is obtained for each field configuration 0(x). This 
means that S* is a functional of the fields involved, which take the role of 
the paths in the mechanical system. 

The simplest case is a neutral scalar field 0(x) G IR. If this field 
describes free scalar particles of mass m, its Lagrangian read^ 

£(0, d,<P) = -d,^{x)d^<P{x) - —ct>{xf . (1.6) 

Assembling the Lagrangian C only by covariant terms — as it is the case 
in eq. (11.61) — ensures that we are dealing with relativistic field theories. 

In classical field theory the configuration is determined by again en- 
forcing the variational condition 5S = 0. For a neutral scalar field, this 
implies 

dcf) d{d^(j)) 

which translates for the Lagrangian (11. 6p into the Klein-Gordon equation 
of motion for the scalar field, 

[df.d^' -m^](j){x) = . (1.8) 



-'^For convenience, we set the speed of light to c = 1. 
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In simple situations, the variational principle and the boundary condi- 
tions fix the classical field configuration everywhere in space-time. 

As a well-known example, electrodynamics deals with vector fields 
Afj,{x), which represent the electromagnetic potentials. The Lagrangian 

^ '^F^pF^ -\- j^A^ , Fij^p dfj^Ay dyA^ , (1-9) 

is constructed from the (gauge invariant) field strength tensor F^^, and we 
added an external, electrically charged current j^{x). Now the condition 
6S = leads to the inhomogeneous Maxwell equations 

d.F^" = f , 

from which we infer that the current classically obeys the continuity equa- 
tion d^j'^ = 0. (On the other hand, the homogeneous Maxwell equations 
are already encoded in the use of potentials.) 



1.3 Quantum field theory 

The transition from classical field theory to quantum field theory can 
be performed in analogy to the quantisation of the mechanical system 
in Subsection 1.1. Since the role of paths in that case is now taken 
by configurations, we quantise the field theoretic system by including 
contributions of all possible field configurations. To render such a huge 
summation well-defined, we introduce again a discretisation. Since the 
fields take their values in each space-point x, we now need a space-time 
lattice, which we choose to be hypercubic, and we denote the lattice 
spacing again by a. Thus the lattice consists of the sites 

^ G ^ , V/i| . (1.10) 
a J 

If we stay with the example of a neutral scalar field, then all the config- 
urations are summed over as follows, 

/oo n 
d(t>.... (1.11) 

X fc lattice 

In this summation, we are going to attach a phase factor exp(-|S'[0]) to 
each configuration, similar to eq. (11.31) . On the right-hand-side we indi- 
cate again the continuum limit, the details of which will be of prominent 
interest in this work. 

A configuration which corresponds to the lowest possible energy is 
denoted as a vacuum Vt. Similar to eq. (11.91) we add an external source 




10 



field J{x), which now couples to the field Then the vacuum-to- 

vacuum transition amplitude is defined as 

Z[J] = {Q\n)j = j'^'P exp + J0)) , (1.12) 

where we use continuum notation, and Jcf) = J d'^x J{x)(j){x). 

Let us assume for simplicity the solution of the equation 55 = to be 
unique. Then it is again the vicinity of this classical configuration which 
contributes in a dominant way (on an action scale where h is small) ; also 
here the contributions at large \6S\ are mostly washed out by the rapidly 
rotating phase. 

The convergence of the sum over the configurations can be accelerated 
drastically if we perform a Wick rotation t —it to arrive at Euclidean 
space. There we denote a point as x = (x, 0:4), X4 being the Euclidean 
time, and the above quantities turn into 

£e(0,5/.0) = ^df,^df,(f) + V{(f)) , Se[4>] = j d^xCEi(p,d^,(p) , 

Ze = Ze[J = 0] = Jv<f)expl^- ^5e[0]) . (1.13) 

Ce and Se are the Euclidean Lagrangian and action. V{(f)) is some 
potential, which is — for instance — quadratic in the free case, as we 
saw in eq. (11. 6p . In Euclidean space we only write lower indices, and 
doubled indices are summed over from 1 to 4 with the metric tensor 6^^- 

Now the contributions by configurations deviating from the action 
minimum are suppressed exponentially, which speeds up the convergence 
of the functional integral tremendouslyo This property is highly welcome 
if we try to evaluate functional integrals approximately by summing over 
a small but (as far as possible) representative set of random configura- 
tions. This is the method used in numerical simulations, which we will be 
concerned with later. If conclusive simulations are feasible, they provide 
in most cases the only access to functional integral results beyond pertur- 
bative, semi-classical or effective approximations, i.e. to actual functional 
integrals at finite interaction parameters. 

In the terminology of statistical mechanics, Ze is a partition function. 
Then h takes a role analogous to the temperature, which controls the 
extent of field fiuctuations around an action minimumjf] In the limit 
h —>■ only the latter contributes (the system "freezes" to the classical 

^Of course, the Wick rotation could have been performed ear her in quantum me- 
chanics, where it accelerates the convergence of the path integral as well. 

^In this sense, the variation of h does have a realistic interpretation, although it is 
fixed in Nature. 
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configuration), so this limit leads back to the classical field theory of 
Subsection 1.2. Once more there is an analogy to the point mechanics in 
Subsection 1.1. 

In quantum field theory, the fiuctuations around the vacuum are es- 
sential; they record the occurrence of particles, deviating from a vacuum 
state Q. 

In view of the statistical interpretation, we can build expectation 
values, and these are the quantities that contain the physical information. 
The vacuum expectation value of some observable 0{(f)) is given by 

{0{<p)) := {n\o{m) = ^/ ^'^^('^) (^('^t'^])) ' (1-14) 

so that eq. (11.121) fixes the normalisation (1) = 1. In particular, the 
Euclidean 2-point function takes the form 



J=0 



= ^yp00(x)0(y)exp(-i5E[0]) . (1.15) 

Here we assumed the condensate (or 1-point function) (0) to vanish, 
hence G2{x — y) coincides with the connected correlation function (with 
the general form {(j){x)(j){y)) — {(j){x)) {(j){y)) ). It characterises the corre- 
lation over a temporal separation At = — y4^ and a spatial distance 
X — y. If one Fourier transforms the distance, one usually obtains an 
exponential decay in At, 

G'2(p,At) oce-^(^^)^* , (1.16) 

where E{p) is the energy of the particle involved; in particular E{0) is 
its mass. 

Similarly we may extract further information of physical interest by 
evaluating higher n-point functions 

G„(a;(i),...,a;(")) = (0(x«)---0(x("))), (1.17) 

or their connected part, which is often of primary interest, 

. (1.18) 



J=0 



Here all the x^^^ are Euclidean space-time points. 

In the further Sections we will stay in Euclidean space (unless it is 
specified otherwise), and we will from now on suppress the subscript "E". 
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The use of the Euchdean signature is justified because the expectation 
values — which provide the physical observables — can be carried over 
to Minkowski space, if four conditions are fulfilled. These conditions are 
known as the Osterwalder-Schrader axioms [2j. Two of them ("analytic- 
ity" and "regularity") are rather technical, while "0(4) invariance" and 
"reflection positivity" have a physical interpretation. Note also that n- 
point functions in Minkowski space require a time ordering. If we deal 
with Euclidean lattices, we assume first a continuum limit to be taken, 
and then the transition to Minkowski space to be justified. 

From now we use on natural units, h = c = 1, and — when it is 
specified — also lattice units, which set in addition the lattice spacing 
a = 1. In Sections 7 and 8 we identify the spacing in lattice QCD with 
a physical scale, which then attaches physical units to all dimensional 
quantities involved. 

Derivations and details of the basic features that we have sketched 
in this Introduction can be found at numerous places in the literature. 
Due to their established status, we hardly indicated references so far. At 
this point we would like to attract attention to Ref. [3], which covers the 
subjects hinted at in Section 1 with great precision. This also includes an 
explanation of the Osterwalder-Schrader axioms and a comprehensive list 
of references on the functional integral formulation of quantum physics. 

2 Renormalisation Group Transformations 
and Perfect Lattice Actions 

2.1 Block variable transformations 

In Section 1 we have introduced the partition function Z and its func- 
tional derivatives as the quantities of interest. They are well-defined on 
a Euclidean lattice, which restricts the momenta to the Brillouin zone 



i.e. it naturally introduces a momentum cutoff A = vr/a. 

For a neutral scalar field, G IR , the source-free partition functions 
takes the form 



Also the action S[(j)] is affected by the discretisation. The standard form 




(2.1) 




(2.2) 



X 
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for a free lattice scalar field reads 



5Z f 2a2 5Z 



,2 



' C 1 

where fi is the unit vector in yU-direction, and 



a y 
2^) 



2 . ap^ 
— sm — - 
a 2 



P 



d 
^Ji=l 



(2.3) 



(2.4) 



As this modified momentum shows, the lattice structure introduces arti- 
facts on a scale fixed by the cutoff A. 

Now we would like to alter the lattice action in a way that moves the 
cutoff effects to higher energies. This can be achieved by a renormali- 
sation group transformation (RGT) [3] to a new lattice field 0' living on 
a coarser lattice {x'}, for instance with spacing 2a. We can choose the 
sites x' as the centres of disjoint unit hypercubes on the fine lattice {x}. 
Then the action S' for the lattice field ch' can be formulated as 



n 



exp 



-S[<p] - a ^ (^0;, -^^(px 



Xdx' 



(2.5) 



where the sum x G x' runs over the 2'^ sites on the fine lattice in the unit 
hypercube with centre x' . a > and h are RGT parameters, which will 
be commented on below. 

The RGT (12. 5 p leaves the partition function invariant (up to a con- 
stant factor) j§ 



Z' 



n 



-5M 



exp 



a 



xdx' 



Z ■ const. 



(2.6) 



Also the n-point functions are transferred to the coarse lattice without 
any damage, for instance 



0L' (t)'y' ) = {(yY(t)^(yY(t)y 



(2.7) 



xGx' 



y&y' 



^In this case, the artifacts occur in 0{h?), which is generic for bosonic systems. 
^This constant factor does not have any impact on physical properties, since it 



drops out in aU expectation values. 
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If we now consider the situation in terms of lattice units, we set on 
the fine lattice a = 1 and on the coarse lattice a' = 2a = 1. In the 
transition from the former to the latter lattice units, the fields and the 
parameters are re-scaled according to their dimensions. But the use of the 
blocked actions S'[(j)'] guarantees that the lattice artifacts — in particular 
the discretisation errors in the kinetic term — are still those of the fine 
lattice. Hence their scale is A = 2A', in contrast to the standard action 
on the coarse lattice. 

This blocking variable RGT can be iterated, and — for the RGT 
parameter b = 2'^^^~^ — the lattice action converges to a finite fixed 
point \SBS% 

S — >S' — > S" — > ... — >S* . (2.8) 

The scale for the lattice artifacts is unchanged, hence it diverges in lattice 
units, 

A — >2A' — > 4A" — > ... — >oo . (2.9) 
Therefore, S* is free of any cutoff artifacts; it is a perfect lattice action. 



2.2 Blocking from the continuum 

Let us now generalise the blocking factor to ri G {2,3, ... }, so that 
a' = n ■ a. The limit n — >■ cxo means that we perform a blocking from 
the continuum] on the scale of the blocked lattice units, the initial lattice 
appears continuous. Thus we arrive at the perfect action in one single 
step, S ^ S*. The corresponding transformation reads 




(2.10) 



where now (j) is the final lattice field on the sites x, while (f is a continuum 
field, S is the continuum action, and is the unit hypercube (in final 
lattice units) with centre x. Since the RGT does not alter physical 
properties, we see here directly that S* captures continuum physics on 
the lattice, without any possibility for lattice artifacts to sneak in. 

Let us now look at the explicit form of the free scalar propagator 
for the standard lattice action (12.31) and for the perfect action S*[(j)] in 

^For some field of dimension [Mass]''' and a blocking factor n, the corresponding 
parameter multiplying J2xex' ■ ■ ■ for the sake of convergence 

under RGT iterations. This factor compensates the re-scaling at the end of each step. 

^Also the continuum field (p is expressed in the (upcoming) lattice units. 
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momentum space E] (in lattice units) , 



11 1 

:G2(p) standard = " ('/'("P) standard 



1^ / N 1/^. N\ n2(p + 27r/) 1 

2G2(p)perfect = - perfect = 2^ 2vr/)2 + ^ a ' 

where nfp) := TT ^ . (2.11) 

We see that the perfect propagator consists of the continuum propagator 
with an analytic factor 11^ and all 27r periodic copies, plus a constant 
term. The latter vanishes in the limit to a (5-function RGT, o; — > oo. 

The function Il{p) ensures that the sum over the integers converges. 
In the exponential transformation term in eq. fl2.10p we have used a 
step function shape for the integration of ip (1 inside Cx, otherwise). 
This shape could be varied, which implies different forms of the function 
n (without danger for the converges of the sum over the copies of the 
Brillouin zone). For instance, the generalisation to B-spline blocking 
functions is discussed in Ref. [6]. 

The function 11 (p) in eq. (12. lip is also affected by the hypercubic 
structure of the lattice; as an example, we could stay with the step func- 
tion shape and consider a 2d triangular lattice instead, which leads to a 
different function in the enumerator ^ (in this case, the lattice cells to be 
integrated over are the hexagons of the dual lattice). However, it always 
has to be analytic]^ hence it does not affect the dispersion relation to be 
extracted from the perfect propagator The latter coincides indeed with 
the continuum dispersion, which also means that it displays exact and 
continuous rotation invariance (which turns into exact Lorentz invariance 
in Minkowski space). We emphasise that this symmetry can be imple- 
mented in the physical observables (not in the form of the propagator), 
irrespective of the lattice structure. 

In coordinate space we write the perfect action in the form of a dis- 
crete convolution 

^*[0]=^0.p(a;-y)0, , (2.12) 
x,y 

where p{x — y) is the inverse Fourier transform of G2{p)~^y:iecf '^^^ decay 
of \p{x — y)\ is exponential in \x — y\ for any mass m (for increasing 
the decay is accelerated). Explicit examples are shown in Ref. [B]. 
This means that the perfect lattice action is local. Generally, locality 



*We assume the analytic continuation to be inserted at the removable singularities. 
^Of course, one always considers the branch with the lowest energy. 
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is reputed as a vital requirement to ensure that a lattice action has a 
sensible continuum limit. 

However, in d > 2 the couplings in p extend to infinite distances 
\x — in contrast to the ultralocal standard action]^ For practical 
purposes this set of couplings has to be truncated. To this end, we first 
identified the value of a which optimises the level of locality, i.e. the 
rapidity of the exponential decay. This optimal value depends on m; a 
good approximation is [6] 

sinh m — m 

"optimal ("^) ^ ^ , (2.13) 

which is derived from the property that it only couples nearest neighbours 
in d = 1. Then we may truncate p to a unit hypercube — i.e. we 
enforce that couplings p{x — y) only occur if — Vfil < ct, V/x — 
by imposing periodic boundary conditions over 3 lattice spacings. This 
method yields a perfect action in a lattice volume 3^, which we then use 
as an truncated approximation to the perfect action on larger lattices. 
Unlike other truncation schemes, this one guarantees for instance the 
correct normalisation of the couplings. 

Figure [1] illustrates as examples the dispersion relation E{p) (for mo- 
menta p oc (1, 1, 0)) for the free scalar particle of mass m = 2, as well as 
the thermodynamic ratio P//i^ at m = (where P is the pressure and n 
the chemical potential) 1"1 



2.3 Classically perfect actions 

For interacting theories, the perfect action can in general not be com- 
puted explicitly, since this requires carrying out a functional integral. 
P. Hasenfratz and F. Niedermayer ^ suggested a feasible simplification, 
which evaluates the RGT steps in the classical approximation. This idea 
has revived and boosted the RGT method in lattice field theory. In our 
case, this classical RGT step takes the form 

m = T {s[<P]+aj2 (0;, - ^^E'^^)'} • (2-14) 

x' xdx' 

Iteration leads also here to a fixed point — the classically perfect action. 
With a suitable parameterisation ansatz, it can be determined numeri- 
cally to some approximation by inserting a set of configurations 0' on the 

^°Ultralocality means that the couphngs drop to zero beyond a finite number of 
lattice spacings. 

"'^"'^The inclusion of a chemical potential in a perfect lattice action will be commented 
on later in the fermionic context (Subsection 4.1). 
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Figure 1: On the left: The dispersion relation for a free scalar particle 
of mass m = 2 for the perfect lattice action (which coincides with the 
continuum dispersion), the hypercubic, truncated perfect action and the 
standard action. For increasing momenta, the magnitude of the lattice 
artifact rises rapidly for the standard action, whereas it remains modest 
for the truncated perfect action. 

On the right: the corresponding scaling test for massless lattice scalars 
with respect to the ratio between the pressure P and the fourth power 
of the chemical potential fi. As /i increases, lattice artifacts cause a 
deviation from the continuum value = l/(487r^). This deviation 

is large for the standard action, but harmless for the truncated perfect 
hypercube action. 

In both plots all quantities are given in lattice units. 



coarse lattice and performing the minimisation. The parameters which 
are used in the ansatz for the action are then tuned until one obtains 
optimal approximate invariance under this transformation, i.e. an ap- 
proximate classical fixed point action. This procedure is particularly 
promising for asymptotically free theories: there the weak field expan- 
sion of the classical fixed point equation corresponds in leading order 
the (soluble) case of free fields. The pioneering work [7J for this method 
evaluated and simulated a classically perfect action with a large number 
of parameters in the 2d 0(3) model (a non-linear a-model). A subtle 
scaling test (suggested in Ref. [8]) revealed practically no lattice artifacts 
at all down to (^/a ~ 5 (where ^ is the correlation length). This is in 
contrast to the standard action, where scaling artifacts are visible even 
at ^/a ~ 15 [2]. 

Later applications of classically perfect actions include topological 
aspects of the 2d 0(3) model 0, the 2d 0P(3) model [ID], pure SU{2) 
[TT] and SU{3) \T2\ gauge theory in = 4, the two-flavour Schwinger 
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model [13! and finally QCD [H]. 

In Figure [2] we show a comparison — involving classically perfect 
actions — for the scaling of the thermodynamic ratio P/T^ (where T is 
the temperature) for free scalars [6] (on top), and for the static quark- 
antiquark potential fT5] (below). 
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Figure 2: On top: the scaling ratio (pressure) /(temperature) for finite 
numbers Nt of lattice sites in the temporal direction. A decreasing num- 
ber Nt corresponds to a coarser lattice, which amplifies the artifacts, in 
particular for the standard action. (The continuum value is given by the 
Stefan- Boltzmann law, P/T^ = tt'^/90.) 

Below: The (re-scaled) static quark- antiquark potential V{r) at differ- 
ent distances. Wilson's standard formulation is only defined at discrete 
distances and exhibits significant artifacts at short r. The classically 
perfect potential captures all distances and suffers much less from lattice 
artifacts. 

In Ref. [TB] we studied a free scalar particle on a circle (a quantum 
rotor) with a discrete Euclidean time and periodic boundary conditions 
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over a period T. We considered the scaling of the ratio between the 
first two energy gaps and of the topological susceptibility (scaled by the 
correlation length 



E2 — Eq 
El — Eq 



and 



Xt 



(2.15) 



The latter is based on the expectation value of the squared winding 
number u, which is the simplest case of a topological charge. These 
results are plotted in Figure O 
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Figure 3: The scaling behaviour of a free scalar particle on a circle. We 
show two scaling quantities as functions of the correlation length in lattice 
units: the ratio between the first two energy gaps (on the left) and the 
topological susceptibility (on the right). 

It is remarkable that also the continuum topology can be represented 
exactly on the lattice, thanks to the formulation with perfect actions 
and operators. Generally, we build (classically) perfect operators from 
the lattice fields obtained by (classical) blocking [15]. 

In contrast, for the standard action it is not even obvious how to 
define topological sectors (since all lattice configurations can be contin- 
uously deformed into one another). We use it here with the geometrical 
definition of the topological charge [17], which is the best option, but 
we observe strong scaling artifactso On the other hand, the perfect 
formulation keeps track of each detail in the intervals between the dis- 
crete time points, since it emerges from blocking transformations. This 

^^In Figure [3] we use the continuum correlation length as the scale. If one inserts 
instead for the standard action the correlation length as obtained from standard action 
simulations, the artifacts are reduced, but the hierarchy in the scaling quality persists 

m- 
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means that any winding number between nearest neighbour time sites 
is included as a possibihty in the expectation value (z/^). (Of course, a 
large number of windings is strongly suppressed by the kinetic term in 
the exponent of the Boltzmann factor, cf. eq. (11.131) ). The classically 
perfect action still approximates the continuum value of Xti to a very 
good approximation. 



3 Fermions 



3.1 The Dirac equation 

For convenience we temporarily return to Minkowski space for the Sub- 
sections 3.1 to 3.3, which contain introductory remarks on fermions. 

Let us go back to quantum mechanics as a renewed starting point. 
Taking the relativistic energy-momentum relation E'^ = p'^ + m? as a 
guide-line, one arrives at an obvious ansatz for a relativistic Schrodinger 
equation, 

[d^,d^ - m^] ^ = . (3.1) 

This is the Klein-Gordon equation, which we already encountered in eq. 
f ll.8l) in the context of classical field theory. An apparent problem with 
it, which worried the pioneers of quantum mechanics, is the occurrence 
of negative energies. P.A.M. Dirac wanted to avoid them by linearising 
this equation with the ansatz 



[i-ff.d^' - m] ^ = 



(3.2) 



In order to reproduce the relativistic energy, the coefficients 7^ have to 
obey the anti-commutation relation 



= '^9t,u , 5- = diag(l, -1, -1, -1) 



(3.3) 



Therefore, these coefficients 7^ in the Dirac equation (13.21) have to be 
(at least) 4x4 complex matrices in ci = 4. Thus the spinor has four 
components!^ 

Actually this linearisation does not overcome the negative energy eigen- 
values. Nevertheless this ansatz was extremely successful; for instance. 



(3.4) 



13 



In two dimensions, we can live with 2x2 matrices 7^ and 2-componcnt spinors. 
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it led to the prediction of the positron just before its experimental dis- 
covery in 1931. In fact, the spinor \1/ captures a spin-1/2 particle plus its 
antiparticle. 

Later on, relativistic quantum mechanics considered the Dirac equa- 
tion appropriate for fermions, and the Klein-Gordon equation for bosons. 



3.2 Fermionic field theory 

In the functional integral formulation of fermionic field theory, the Dirac 
operator is still present as the central ingredient in the Lagrangian. For 
free fermions of mass m, the partition function and the action are written 

as 



Si^,"^] = J d^x^{x)[i^^d>' -m]^{x) , (3.5) 

where ^ (a;) = (^/'i(a;), ■?/'2(x), ^/'3(x), '?/'4 (a;)) and \I'(x) are spinor fields. Ap- 
plication of the variational principle 6S = leads to the Dirac equation 
(13.21) for and to the adjoint Dirac equation 

^ [iji!d + m] = . (3.6) 

In the light of the Spin-Statistics Theorem, fermion field components 
anti-commute, hence one describes them by Grassmann variables. A 
set of Grassmann variables {?7i}, i = 1, 2, . . . (as it is used here for the 
components of ^ and of \1/ in a specific point x) obeys the relations 



{'Hu '7j} = , j dr]i r]j = 6ij . (3.7) 

A striking difference from the Dirac algebra (13. 3p is the property rjf = 
0. The integration rule is motivated by the analogy to the translation 
invariance of the real, unbounded integral. The Grassmann integral has 
no bounds, and its effect is equivalent to differentiation. It provides the 
basis for the functional integral in eq. (13.51) jl9] . which we will make 
explicit in Subsection 3.4. 

Interactions can be included for instance by adding a 4-Fermi term 
(^'(x) \l/(x))^ to the Lagrangian, which we will consider in Subsections 
3.5, 4.3 and 5.2. Another type of interaction is generated by coupling the 
fermions to a gauge field Afj, through a covariant derivative, which turns 
the Dirac operator and the partition function into 

D{A) = t^^[d''-gA''{x)]-m, (3.8) 
V^V^VAexp (^i j d^x^{x)D{A)'^{x) + iS[A] 
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S[A] represents the pure gauge action; it could be for instance the Abehan 
gauge action obtained by integrating the Lagrangian fll.9l) . In that frame- 
work, the term g^'jfj,'^ takes the role of the external, charged current 
(and g is the gauge coupling). Fermionic n-point functions are defined in 
analogy to the bosonic case (see Subsection 1.3), but the order matters, 
of course. 

3.3 Chiral symmetry 

Due to the anti-commutation rule (13. 3p . the matrix 

75 := ^7o7i7273 obeys {75, 7^} = and 75 = 1 • (3-9) 
Therefore, the operators 

P±:=l(]l±75) (3.10) 

are complementary projectors (P| = P± , P+ + P- = 1). They can be 
used to decompose the spinor fields into their so-called left-handed and 
right-handed parts, 

^L,i?(x) = P±^(x) , ^lA^) = ^ix)P^ . (3.11) 

In these terms, the fermionic part of the Lagrangian in eq. (13.81) reads 

C = ^LDiAU^o'^L + ^RD{A)m=o'^R - mi^L'^R + ^r^l) . (3.12) 

In the chiral limit m ^ the left-handed and the right-handed parts 
decouple completely. This property corresponds to the relation 

{^™=o,75} = 0, (3.13) 

which manifests itself in a global symmetry, namely the invariance of C 
under the "chiral rotation" 

i|r ^ ^ e^"^5 , ^ ^ e^"^5 ^ (3.14) 

for an arbitrary parameter a. 

Obviously, the term that enters the Lagrangian (13.12^ for a fermion 
mass m 7^ breaks the chiral symmetry explicitly; the chiral rotation 
f l3.14p transforms the mass term as 

m^x^f ^ m^^e^'"^^^ . (3.15) 
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In general, global symmetries — such as the chiral invar iance — are 
only realised approximately in NatureH hence a breaking by a mass 
term is not necessarily a problem. However, by the introduction of gauge 
fields one arrives at local symmetries, and they have got to be exact. In a 
vector theory, the gauge fields couple in the same way to the left-handed 
and to the right-handed fermions. This is the case for the gluon fields in 
QCD. Then the fermion mass is allowed — quark masses can be inserted 
into the QCD Lagrangian. 

On the other hand, the electroweak sector of the Standard Model is 
an example for a chiral gauge theory, where the gauge fields couple in 
different ways to the left-handed and to the right-handed fermions. Then 
we have to require their invariance under independent transformations, 
which forbids explicit mass terms in £. In that framework, the masses 
of fermions (and also those of gauge fields) can only be generated dy- 
namically. It takes Yukawa couplings to a Higgs field and spontaneous 
symmetry breaking to arrive at massive quarks and leptons (and massive 
gauge bosons and Z°)0 

3.4 Fermions on a Euclidean lattice 

We return to Euclidean space, where the 7-matrices obey 

{7,i,7u} = 25^^ , 75 := 71727374 , {7m'75} = , 75 = 1 • (3-16) 

We choose them to be Hermitian. We write a bilinear fermionic lattice 
action — such as the action for free fermions, or for fermions interacting 
through gauge fields — in the form 

N 

S =Y^ ^^Mij^J = . (3.17) 

Here the components run over all the lattice sites, and on each site 

over all internal degrees of freedom (spinor and fiavour indices, and for 
instance in QCD also colour indices). It is easy to see that the partition 
function is given by the celebrated fermion determinant, 

r ^ - 

^ i=i 

^''An exception is CPT invariance, which is assumed to be exact [20] (this was first 
postulated by W. Pauh in 1955). Its breaking would imply the violation of Lorentz 
invariance |21j . which appears as a global/local symmetry in special/general relativity. 

^'"^We do not consider renormalisation effects at this point. 
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This expression also attaches an exphcit meaning to the Grassmann func- 
tional integrals]^ (The order of the single Grassmann integrals matters, 
since the rule in eq. (13 .Tp refers particularly to the innermost integral.) 

As an example, we consider a free, massless fermion with the Eu- 
clidean continuum action 

= J d^xi/j{xh^d^^/j{x) . (3.19) 

On a lattice with unit spacing (a = 1) the simplest discretisation ansatz 
reads 

S[^,^] = ^^r^,7^.1(vl/,,^^,-vI/,._.) 

X 

= j^d^p^i-p)t7^,smp^<il{p) . (3.20) 

This formulation is known as the naive lattice fermion. As the name 
suggests, there is a serious problem with it: the propagator 

G'2(p)„aive = -. ^ (3.21) 

has inside the (first) Brillouin zone (12. ip not only the physical pole at 
p = 0, but it has poles whenever p^ e {0, tt}, = 1 . . . 4. Hence there are 
16 poles (in general, 2°^ poles) instead of the one that we have ordered. 
This effect is known as the fermion doubling problem — it is due to the 
occurrence of a linear derivative. In fact, these doublers distort physical 
properties regardless how fine the lattice might be, hence this formulation 
cannot be applied. For instance, among the 16 species the chiralities 
are equally distributed, which makes it impossible to construct a chiral 
gauge theory [22] • Moreover, the trouble also affects vector theories, 
since these species contribute to the axial anomaly with alternating signs, 
hence doubled lattice fermions cannot reproduce a non-vanishing axial 
anomaly either |23j . 

So we have to consider further options for the lattice action 

S = J2^-D^y'^y ■ (3-22) 

We recall that locality is in general a requirement for a controlled contin- 
uum limit (in view of the extension to the interacting case). In coordinate 
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It is entertaining to compare this result to the expression for a complex scalar 
field, J YlidRe^idlm^i exp(-$^M$) cx (detA/)^^, or l/VdetM for a neutral scalar 
field. 
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space, a local lattice Dirac operator D has to be bounded ai 



\D^y\ < cie-'=2l^-2'l , ci, C2 > . (3.23) 

In momentum space this means that D{p) = G2{p)~^ has to be analytic. 

It is not easy to find a satisfactory solution to the doubling problem. 
This statement was made precise by the Nielsen-Ninomiya No-Go The- 
orem [22] • Putting aside technical details^ it essentially states that an 
undoubled lattice fermion cannot be chiral and local at the same time. 

For an intuitive and simplified illustration, we write a rather general 
ansatz for a chiral lattice Dirac operator for free fermions, 

D{p) = tp.iphf, , p^ip) =P^ + 0(p2) . (3.24) 

The leading momentum order of p^ is required by the correct continuum 
limit (which is determined by small momenta in lattice units). We may 
consider the specific momenta p = (pi, 0, 0, 0), so that 

D(pi,0,0,0) = «pi(pi)7i , 

with the physical zero at pi = 0. Since 2tt periodicity is mandatory, and 
since locality requires an analytic function pi(pi) = Pi + 0{pl), at least 
one additional zero (generally: an odd number of them) is inevitable 
inside the Brillouin range pi G (— 7r,7r], even if one deviates from the 
naive form pi(pi) = sinpi. 

Many suggestions have been made to circumvent this problem by 
breaking one of the desired properties on the lattice, hoping this would 
not affect the continuum limit. We do not review all these efforts, but 
we mention as an example the SLAC fermion [23]. In the above conside- 
ration, it sets pi(pi) = pi in pi G (— 7r,7r], which is then periodically 
continued (and the same for the other momentum components). Due to 
the jumps at the edges of the Brillouin zone this formulation is non-local. 
The hope to get away with this was crushed by Karsten and Smit, who 
showed that this formulation is inconsistent at the one-loop level of gauge 
theory, where it fails to reproduce Lorentz symmetry in the continuum 
limit [2S]@ 



^^Different definitions of locality appear in the lattice literature, but the condition 
of an exponential decay — which we referred to already for scalar fields — is the 
relevant one, because it guarantees a safe continuum limit. 

""^^The proof requires some additional assumptions — like lattice translation invari- 
ance — but they are not especially tricky. 

^^However, this conceptual problem at the one-loop level is specific to gauge inter- 
actions. The SLAC fermion may still be in business for instance in supersymmetric 
spin models without gauge fields [26j . 
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The standard lattice fermion formulation, which has been used most 
in simulations — in QCD in particular — was put forward by K.G. 
Wilson in 1979 [27]. The free Wilson operator reads 



,xy ~ 



'x,y+jl 



(3.25) 



Wilson added the term in the second round bracket to the naive form 
that we considered before. This term represents a Laplacian operator, 
which is discretised in the simplest wayo In fact it avoids the fermion 
doubling by sending the doublers to the cutoff energy. There is no doubt 
that this operator is local, and the Wilson term is 0(a) suppressed, so 
one could hope that it does not distort the continuum limit. 

However, due to this extra term the chiral symmetry is broken expli- 
citly, {-Dw,m.=o 5 75} 7^ 0. As interactions are switched on, this causes 
numerous problems. In particular, a gauge field can be added as a set of 
link variables in the gauge group, which provides invariance under gauge 
transformation of the matter fields on the sites. One often writes this 
compact link variable as 



which indicates a connection to the (non-compact) continuum gauge field 
A^. For non-Abelian gauge groups this exponential is formulated as a 
path ordered product [2H]. Such a gauge field suppresses the terms 5x,y±^j, 
in the Wilson term, but not its last entry 25^ y. This different treatment 
gives rise to additive mass renormalisation. If one tries to approach the 
chiral limit, where the renormalised fermion mass vanishes, one has to 
fine tune the bare mass to some value m < 0, which compensates for the 
additive renormalisation. 

A further (related) inconvenience for interacting Wilson fermions is 
that the lattice artifacts can appear in 0(a) alread}0 (unless one adds 
another term to cancel the 0(a) artifacts — following Symanzik's pro- 
gram — which requires fine tuning again [2U]). 

Another formulation, which has been considered standard over the 
past decades, and which is regularly applied in simulations, is known as 
the staggered fermions (or Kogut-Susskind fermions) [30]. An elegant 

^°The Wilson term can also be multiplied by some independent coefScient (the 
Wilson parameter), but this generalisation is not particularly fruitful. 
^-'^For the free Wilson fermions, the scaling artifacts are of O(a^). 



U^,,x = exp 



i I dy^A^iy)) e { gauge group } , (3.26) 
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way to construct them starts from the naive action on a unit lattice, 
1 

X ^=1 

(3.27) 

and performs the substitutions [31] 

= ^.7r • • • 7? , K = 1? • • • 7r^x . (3.28) 

This leaves the mass term invariant and renders also the kinetic term 
diagonal in the spinor space. Hence one may reduce the transformed 
spinors to a single component x, and one obtains 

1 

S% X] = [2 5Z ^/^'^ {XxU y,^xXx+{, - Xx+fiUlccXx) + mXxXx , 

X M=l 

r^^.x := (_l)^i+^2+-+^M-i . (3.29) 

This structure distinguishes two sublattices by the criterion if Yl'^=i 
is even or odd, i.e. by the sign term 

e{x) = (_i)-i+-+-<* . (3.30) 

The link variables Ufj,^x always connect sites belonging to distinct sublat- 
tices. For 171 = 0, the action (13.291) is invariant under the transformations 

Xx^e'^^^hx, Xx^e'^^^hx, (3.31) 

which amounts to a remnant chiral symmetry U{l)e ® ^(l)o; where e{x) 
adopts the role of 75 (the subscripts refer to the even/odd sublattice). 
The 2"^ components on the corners of the disjoint unit hypercubes are 
nowadays denotes as "tastes" (the earher literature also called them 
"pseudoflavours" ) . As long as one finally assembles exactly 4 flavours 
from them, this remnant symmetry is sufficient to avoid additive mass 
renormalisation and 0{a) scaling artifacts. 

Recently it became fashion to try to build single flavours with stag- 
gered fermions by taking the fourth root of the fermion determinant 
fl3.18p . However, it is likely that this formulation is non-local, and — 
even if someone is willing to accept that — additive mass renormalisa- 
tion sets in again (see e.g. Refs. [32]). The question if this formulation 
might provide correct results even if it is non-local is still debated [33] . 
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3.5 Are light fermions natural ? 

We would like to stress that keeping track of the chiral symmetry is not a 
specific problem of the lattice regularisation. It should rather be viewed 
as a generic and deep problem, which plagues other regularisations as 
well. For instance, in dimensional regularisation [34j one performs com- 
putations in 4 + e dimensions (in the sense of distribution theory) and 
sends e ^ at the end. This is the most popular regularisation scheme 
for perturbative calculations, but it is a longstanding problem to find a 
generally rehable rule for handling 75 = Y['l=i 7/^ '^^ regularised level 
(or ^n)!=o7M Minkowski signature), generalising eq. (I3.16P (resp. 

(13.91) ). A careful analysis of this issue can be found in Ref. [35j . 

On the conceptual level, this observation means that the existence 
of light fermions in our world actually appears to be unnatural. Nature 
must be non-perturbative, so we can only refer to a non-perturbative 
regularisation scheme when addressing this question, which means es- 
sentially the lattice@ It is possible to formulate light lattice fermions 
— for instance light quarks in lattice QCD — but only with tedious 
and sophisticated constructions (see Section 7), which do not appear to 
mimic a conceivable mechanism in Nature. What could be an acceptable 
mimic is something very simple like the Wilson fermion (13.251) . which, 
however, pushes the fermion mass to the cutoff scale due to a strong ad- 
ditive mass renormalisation — unless a negative bare mass is fine tuned, 
which appears unnatural again@ 

Nevertheless there do exist in particular two light quark flavours, 

rriu, rrid < Aqcd • (3.32) 

The question how this is realised in Nature is a hierarchy problem that 
is not properly understood. The situation is different in pure Yang-Mills 
theory, for example, where (regularised) glueball masses can be made 
arbitrarily light thanks to asymptotic freedom and the absence of additive 
mass renormalisation. But for quarks it does not work in this simple way, 
due to the problems to keep track of an approximate chiral symmetry in 
a regularised system, such that the fermion mass is far below the cutoff. 
In a broader framework, this hierarchy problem raises the question why 

^^A conceivable alternative might be the formulation on a "fuzzy sphere" [55] . 
However, even simulations of models without fermions 37J show that it is not obvious 
to recover the desired continuum limit in these formulations. 

^•^In the light of the properties mentioned in the last paragraph of Subsection 3.4, the 
staggered fermion formulation cannot really be considered a solution to this problem 
either. 
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hadron masses are far below the Planck scale, and therefore why we do 
not just consist of gluons. 

In Ref. |38] we studied the question if this problem could be solved 
(qualitatively) in a brane world. As a toy model, our target theory was 
the 2d Gross-Neveu model [39] with the (continuum) action 



d X 



(3.33) 



where we suppress the flavour index 1 . . . A^. It enjoys a discrete chiral 
Z{2) symmetry 



(^L, ^l) ^ ±(^L, ^^l) , i'i'R, ^r) ^ T(^h, ^r) 



(3.34) 



(the chiral components are defined in eq. (I3.1ip ). With an auxiliary 
scalar field $ the action fl3.33p is equivalent to 



^, $] = / d^x h-f^,^^'^J - + 

L Zg 



(3.35) 



as we see by integrating out the field The sign of $ flips under a Z{2) 
chiral transformation (13.341) . 

In the limit oo, $ freezes to a constant [39j, and ^, \1/ can 

be integrated out. The resulting fermion determinant gives rise to an 
effective potential. 



(3.36) 



In a large volume V, the minima ±$o of Vcs obey the gap equation 



1 1 



dk 



k 



{k = Jkl + kl) . 



(3.37) 



At weak coupling <^ 1 we are dealing with a cutoff A2 ^ $0 and 

m = $0 = A2 e"^/^ (3.38) 

represents the fermion mass, which is generated by the spontaneous 
breaking of the Z{2) symmetry (13.341) . The exponent in eq. (I3.38P ex- 
presses asymptotic freedom. 



Let us proceed to three dimensions, where the action 
S[^,^] = ld'x[^j^d^^ + ^jsd3'^-^{mY] , (/x = l,2) (3.39) 



30 



still has a Z{2) symmetry, 

(^^^,^^)|(,^,,3) ^ t(^^,^^)|(,^,_,3) , (3.40) 

which turns into the discrete chiral symmetry (13.341) after dimensional 
reduction. The 3d gap equation reads 

1 1 



and for a cutoff A3 ^ $0 one identifies a critical coupling Gc = n'^/A^. 
At G > Gc we are in a phase of broken Z{2) symmetry with 

whereas weak coupling (G < Gc) corresponds to a symmetric phase 
($0 = 0). 

Canonical dimensional reduction from 3 to 2 dimensions works in the 
usual way if we start from the 3d symmetric phase, and it leads to light 
2d fermions |38]. However, this is not satisfactory in view of our mo- 
tivation: for instance a non-perturbative treatment at finite (on the 
lattice) should not start from the symmetric phase, because this just 
shifts the problem of fine tuning to = 3. Therefore we focus on dimen- 
sional reduction from the broken phase. 

We denote the (periodicity) extent of the third direction by /3, and 
^ = 1/m is the correlation length. Starting from the 3d broken phase, the 
limit ^^™Q P/^ = 21n(l + -\/2) does not provide light fermions. Hence we 
proceed differently and generate a light 2d fermion as the ^3 = mode 
on a brane. For the latter we make the ansatz $(^3) = $0 ^§^($02^3)5 
which is inspired by Refs. |10]. We choose X2 as the time direction, hence 
the Hamiltonian reads 

H = 72[7i9i + 73^3 - Hxs)] . (3.43) 

The ansatz \l/(x3)e*'^^^^e~*'^* (and the chiral representation for 7j) reveals 
one localised eigenstate of H, 

with energy £"0 = —ki > 0, i.e. a left-mover. (On an anti-brane — $(0:3) 
one obtains a right-mover with Eq = ki > and exchanged components 
in ^'o(a;3))- 
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In addition there are bulk states (not localised in X3), 



with E = ± Y /c^ + $Q , which form together with an orthonormal 
basis for the 1-particle Hilbert space. 

To verify the consistency of the brane profile we have to consider the 
chiral condensate — \l/o does not contribute to it, and if we sum 
up the bulk contributions of ii^ < we reproduce exactly ^{xs), which 
confirms the self-consistency of this single brane world. 

In addition we are free to fill some of the \l/o states. Those with 
Eq < $0 are confined to the (1 + l)-d world, whereas states with Eq > $0 
can escape in the 3-direction. For the low energy observer on the brane 
this event appears as a fermion number violation. 

We now want to include both, and "^n, to be localised on a brane 
and an anti-brane, and (3 now denotes their separation. For the corre- 
sponding profile we make the ansatz 

<I>(X3) = $o(a[tgh+ - tgh_] - 1) , 
tgh± := tgh(a$ok3±/5/2]) , ae [0,1] . (3.46) 

The ansatz for a bound state with the same form as on single branes, 
Mxs) = c f cosh"|(a$o[a:3 - ^/2]) \ 

implies the condition tgh(a$o/3) = a. Hence the parameter a controls 
the brane separation, such that a ^ and a 1 correspond to /3 ^ 
and /5 — > 00, respectively. 

The Dirac equation in this background still has an analytic solution, 
which is given by the ansatz fl3.47p with 



1 / a$ 
c = - 







2 V ^0(^0 + ki) 



Eo = ±^kl + m^ 



ai = -i{Eo + ki) , ^2 = m = Vl - a2 $0 • (3.48) 

The resulting "^qIx^) represents a Dirac fermion with components 

localised on the brane resp. the anti-brane. For a fast motion to 
the left (right) we have < E'o — —ki {+ki), so that the lower (upper) 
component dominates. This situation is sketched in Figure HI The mass 
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m measures the extent of the L, R mixing. The hmit a — does not 
provide a hght fermion (m = $0)5 but the opposite hmit a — > 1 achieves 
this, since the L, R mixing is suppressed as 

m~2$oe-^*o. (3.49) 

Counter-intuitively, large (3 imphes ^ ^ (3 and therefore dimensional 
reduction. A low energy observer in d = 1 + 1 now perceives a point-like 
Dirac fermion composed of L— and -R— modes. On the other hand, a 
high energy observer in (i = 2 + 1 refers to the scale $0 (the 3d fermion 
mass) and observes a Dirac fermion with strongly separated L— and 
i?— constituents. 

Also the bulk states can be determined analytically, 

eifc3X3 / ^ fci) [a$otgh_ - ik^] 



U := 2E{E + k^){h + aHl) , E = ±^kl + kl + ^ 

Summing up again their E <() contributions to —^>^ yields 
G 



N 



dk^^k^^k. 



= <f(x3)+C, (3.51) 



i.e. the desired result up to a term C, which is given explicitly in Ref. 
[38] . It has to be cancelled by occupying bound states in \E'o, which do 
contribute this time. This requires all the bound states with energies 
Eq < Ep to be filled. The Fermi energy turns out to be Ep = ^q, 
i.e. exactly the threshold energy for the escape into the third dimension. 
Hence this brane anti-brane world does contain naturally light fermions, 
but it is completely packed with them, so its physics is blocked by Pauli's 
principle. 



Since this brane anti-brane world is not topologically stable, we also 
checked if the brane and anti-brane repel or attract each other, which 
could lead to disastrous scenarios. However, it turns out that the brane 
tension energy per fermion does not depend on the brane separation, so 
this toy world is indeed stable f38]. 

Finally we studied the possibility of adding a fermion mass term 
M^f\Ef to the Lagrangian, so that the Z{2) symmetry is also explicitly 
broken in ci = 3 (which is actually realistic for a lattice formulation at 
finite A^). This lifts the degeneracy of the minima of V^ff($). If we still in- 
sert the profile fl3.46p . the condition for — ^\E' requires the bound fermion 
states to be filled even beyond $0) hence in this case there is no stable 
configuration at all. 
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One might also start from the symmetric phase and add a mass term 
to construct a somehow natural starting point. However, such a mass 
is simply inherited by the dimensionally reduced model (while the cut- 
off keeps the same magnitude), hence this does not solve the hierarchy 
problem under consideration. 

We drop this mass term again and summarise Subsection 3.5 by re- 
peating that the construction of naturally light fermions is basically suc- 
cessful, but unfortunately this world does not enjoy any flexibility for 
physical processes. However, we assumed translation invariance in the 
2d world so far. That symmetry may be broken at sufficiently large 
chemical potential, so that the chiral condensate prefers a kink anti-kink 
pattern, rather than a constant [Hj. This could possibly provide the 
missing flexibility for a lively brane world of that kind. 



4 Perfect Actions for Lattice Fermions 
4.1 Free fermions 

In the previous Section we have described the severe conceptual difficul- 
ties with the formulation of fermions on the lattice. We now proceed to 
the application of the RGT technique — described in Section 2 — to 
lattice fermions. This is going to reveal how the perfect action handles 
— and solves — the problems of species doubling and chiral symmetry. 

We start with the free fermion and apply immediately the blocking 
from the continuum (introduced in Subsection 2.2), which is most efficient 
for analytic calculations. In analogy to eq. fl2.10p we now relate lattice 
spinor fields to their counterparts in the continuum, 

f d'^uipiu) , "^^r^ 1 d'^v^iv) . (4.1) 

This relation is imposed by the RGT, which leads to the perfect lattice 
action S'['l', \E'] for free lattice fermions. 



exp I - 5^ [fr^ - /" d''ui,\u)\{R-'ri^[^- I d'vi,^{v)]] , {A.2) 

where s['?/','?/'] = / d'^uijj{u)['y^d^ + mlip^u) is the continuum action (cf. 
eq. (13.191) ). On the lattice (of spacing a = 1) we arrive at the following 
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Figure 4: We show the brane anti-brane profile Iji3.46\ ) in the 3-direction 
of our toy brane world JSE]. The 2-direction is the (Euclidean) time, and 
the left-handed (right-handed) fermion moves to the left (right) in the 
one spatial direction inside the brane world. In the plot above it does 
so with a momentum pi = —Am and below with pi = 2m, where m is 
the fermion mass. The latter arises from the communication between the 
left- and right-handed components, which are localised on the brane and 
anti-brane. As they drift apart, m decreases exponentially in the distance, 
hence a low energy observer on the brane perceives them on top of each 
other as a point-like Dirac spinor. 



perfect action and propagator 

s[^,% = j^d'p^{~p)G{prH{p) 



2 



Gip) = Z - "'"VmV (4.3) 
where the function 11 (p) is defined in eq. (12. lip (also here it could be 
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generalised). This formula has been computed in various ways 
[T5]l^ It incorporates the continuum propagator, its periodic copies and 
the blocking term, in full analogy to the perfect action for free scalars in 
eq. dm]). 

Let us now discuss the role of the blocking term, which we have 
generalised from the constant l/a in the scalar case (eq. fl2.5p ) to the 
form K^Jy. For sure we have to require R to be local. Thus it cannot 
disturb the pole structure of G{p). Hence the formulation is free of 
doublers, and the dispersion relatioij^ coincides with the continuum. 

In the limit R ^ we perform a S-function blocking, as we mentioned 
for the scalar fields before. Then the relations (14.1 p turn into equations. 
In this case (or more generally, whenever {i?, 75} vanishes), G{p)m=o — 
and therefore also the Dirac operator D{p)m=o — anti-commutes exactly 
with 75. Then we have chirality, i.e. invariance under the global trans- 
formation (13.141) . just as in the continuum. Hence the question arises 
in which way a contradiction to the (mathematically rigorous) Nielsen- 
Ninomiya Theorem [22] is avoided. The answer is that in this case the 
Dirac operator is non-local [121113]: it does not decay exponentially, but 
only as ^15j 

D(r)„=o oc . (4.4) 
As soon as we proceed to some non-vanishing, local term R, which obeys 

R75}7^0, (4.5) 

locality is restored. However, this obviously leads to 

{D{p)m=0,l5}y^0 , (4.6) 

hence we do not have chirality in the standard form (I3.14p anymore. 

Still, this breaking of the chiral symmetry can only be superficial: 
we know that the RGT does not distort any physical properties, hence 
the chirality of the continuum must be preserved in the physical obser- 
vables, despite the relation (14.60 . as we emphasised at numerous occasions 
O H?! HH]. Therefore this must be a specifically harmless anti- 
commutator. Indeed it gave the crucial clue for a general criterion for 
the form of such a non- vanishing anti-commutator ^S], which is still 
compatible with chiral symmetry in a lattice modified form [19]. This 



^''Later on it turned out that the perfect propagator G{p) was already discussed in 
the Ref. [44]. However, that work was forgotten until it was accidentally re-discovered 
by P. Hasenfratz in 1997 [45], 

25we repeat that one always considers the branch with the lowest energy, cf. Sub- 
section 2.2. 
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criterion is now denoted as the Ginsparg- Wilson relation (since it was 
already mentioned in Ref. [S]), which we will discuss in Section 7. 

We assume R to have the structure of a Dirac scalar, and we move to 
coordinate space. The perfect action for the free fermion is given in the 
form 



Dxy = {G'^)xy = 7mPm(2^ -y) + - y) , (4.7) 



i.e. it consists of a vector term plus a scalar term. We consider the 
local case (14.51) . where Pfj_{x — y) and X{x — y) decay exponentially in the 
distance \x — y\. For practical purposes we need a truncation in these 
couplings, and we follow again the scheme of Section 2: we first optimise 
R, for the case Rl^iy = pS^yS'^^ . An analytic calculation in d = 1 suggests 
the choice 

, , e™" — m — 1 , , 

p(m) = . (4.8) 

Only for this form of p{m) the Id couplings are limited to nearest neigh- 
bour sites, i.e. they take the structure of In c? > 2 couplings over 
all distances are inevitable, but the choice (14.81) still provides practically 
optimal locality, i.e. optimally fast exponential decays of the functions 
p^{x — y) and A(x — y)] this is illustrated in Ref. [15]. 

As a truncation scheme, we computed for this function p(m) the cou- 
plings of perfect actions in a periodic 3^ lattice, and applied these cou- 
plings in larger volumes [50]- This yields the free hypercube fermion (HF), 
which still has the structure of eq. (14. 7p . but now with strictly limited 
supports for the ingredients to the Dirac operator, 

DBF,xy = It^Pfiix -y) + X{x-y) , 

supp[p^(x - y)], supp[A(x - y)] C |x, y \x^, - < 1 , V/ij . (4.9) 

Tables for the explicit couplings for such HFs at various masses are given 
in Ref. [50] ; in Table [1] we display here the HF couplings at m = to an 
extended precision of 16 digits. 

After truncation, the scaling behaviour is still by far superior to the 
Wilson fermion, and also to the so-called D234 fermion which is 
improved to the leading order in the lattice spacing, following Symanzik's 
program. A comparison of the dispersion relations at mass m = and 1 
is shown in Figure [51 We see a striking improvement for the truncated 
perfect action. 

This trend is also confirmed for the thermodynamic quantities plotted 
in Figures and [71 The pressure P at finite temperature T is obtained 
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r 


Plir) 


A(r) 


(0,0,0,0) 





1.852720547165511 


(1,0,0,0) 


0.1368467943177540 


-0.060757866428667176 


(1,1,0,0) 


0.032077284302446526 


-0.030036032105554878 


(1,1,1,0) 


0.011058131255574036 


-0.015967620416694967 


(1,1,1,1) 


0.0047489906005042248 


-0.0084268119917885868 



Table 1: The couplings of the free, massless HF with the parameterisation 
of eq. l!i4-9\ ). p^{r) is anti-symmetric in and symmetric in all other 
components r^, while A(r) is symmetric in all directions. 



by imposing periodic boundary conditions in the time direction over Nt 
lattice points. The corresponding data in Figure [6] are evaluated with 
the formula 



p ^4 ^ 



(27r) 



P4,n=27r(n-l/2)/Aft 



-'^ * n=l 

1 r 1 

/ dpi\ndeiD{p,pi) . (4.10) 

27r J 



In this case we find a deviation from the continuum value P/T' = 
77r^/180 even for the perfect action, because its perfection is designed 
specifically for zero temperature. 

Figure [7] deals with the inclusion of the chemical potential fi. This is 
achieved by the prescription worked out in Refs. [52] . The key observation 
is that starting from any prefect lattice action at /i = and performing 
consistently the substitutions 

^{x,X4) e~^^''^'(f,X4) , ^(f,X4) e''''^-^{x,X4) , (4.11) 

one obtains in fact a perfect action at finite fi. (Also classical perfection 
is preserved under these substitutions.) In our perfect propagator in 
momentum space fl4.3p . this substitution can be implemented by shifting 
P4 — > P4: + ifi- Then one obtains the pressure and the baryon density (one 
third of the fermion density) at T = as 

^(/^) = 777-^1 fin detD (p,p4) -In detD (p,p4 + i/x) 
(27r)* Jb L 

UB = l^Pif^)- (4-12) 

The scaling is then measured by the deviations from the continuum values 
P/n' = 1/(6-7?^) and Ub/ = 2/(97r^). Lattice artifacts are amplified for 
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Figure 5: On top: The dispersion relation for free, massless lattice 
fermions in d = A for spatial momenta p oc (1,1,0) (as an example). 
For the perfect fermion the dispersion coincides with the exact dispersion 
in the continuum, and the HF dispersion follows it closely. The stan- 
dard Wilson fermion deviates strongly at increasing momenta, while the 
Symanzik improved D234 fermion behaves well up to \p\ ~ 1, before it 
hits a doubler coming down from higher energy. 

Below: Dispersion relation for the free HF at mass m = 1. Here we 
show the energy E for various directions of the momentum p (p = \p\) 
to illustrate that they all follow closely the continuum dispersion over a 
sizable part of the Brillouin zone. 



increasing chemical potential jj,. We see that they remain modest over a 
broad range (i.e. up to coarse lattices) for the truncated perfect action, 
in contrast to the Wilson fermion and the D234 fermion. For large 
the fermion density turns into a constant plateau for the usual lattice 
fermion formulations, and one might believe that this is inevitable due 



39 



2.0 



1.5 



1.0 



-0.5 



10 



— exact 

• Wilson 

♦ fixed point 
*D234 

A hypercube 




15 



20 



Figure 6: The ratio between pressure and (temperature) for various 
types of free lattice fermions, compared to the Stefan- Boltzmann law in 
the continuum. The RGT improved actions converge much faster to this 
value for decreasing temperature (increasing Nt) than the Wilson action 
or the D234 action. 



to Pauli's principle. However, the height of this plateau depends on the 
coupling range of the lattice Dirac operator, and it rises to infinity for 
the (untruncated) perfect action — hence the RGT is able to solve this 
problem as well 
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Figure 7: The ratios P/fJ.'^ and tib/ii^, for the pressure P, the baryon 
density ub and the chemical potential n, at zero temperature, for vari- 
ous types of free massless lattice fermions. For the truncated perfect HF 
these ratios converge very fast to the continuum values as fi decreases, in 
contrast to the Wilson fermion and the D234 fermion. 
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4.2 Perfect staggered fermions 

As we pointed out in the previous Subsection, the full standard chiral 
symmetry cannot be preserved in a perfect and local lattice action. On 
the other hand, staggered fermions only have a remnant chiral symmetry 
— see Subsection 3.4 — which raises the question if that symmetry can 
persist in a perfect and local staggered fermion formulation. In fact, 
this U{l)o ® U{l)e symmetry can be preserved under the RGT, if the 
block variables are constructed such that they do not mix any of the 2"^ 
tastes. A corresponding blocking scheme with overlapping blocks was 
first proposed in Ref. [53j. By its iteration we constructed a perfect 
action for free staggered fermions [13j, which does fulfil the [/(l)o®t/(l)e 
symmetry exactly, and which is manifestly local — the Nielsen Ninomiya 
Theorem does not exclude this remnant chiral symmetry^ In d = 2 the 
corresponding free perfect action for the four massless tastes reads 

4 

s[x,x] = (4-13) 

x,y i,j=l 

( a^iv) \ 

-0L2{v) 

a2{p) ai{p) 

\ ai{p) a2{p) / 

where d{p) is a matrix of phase factors, which arrange for the shifts to the 
appropriate lattice sites (it is given explicitly in Ref. [56], which denotes 
it as D(pY^'^). c is an arbitrary (real) RGT parameter, which we tuned 
again for optimal locality. In this case, the analytic optimisation in d = 1 
yields c = 1/2. Ref. [STj discusses the extension of this action to = 4, 
as well as the generalisation to a finite fermion mass m, which fills in 
diagonal elements in the above matrix and changes the locality optimal 
RGT term. 

Also this result can be derived efficiently by blocking from the con- 
tinuum, if the overlapping integration cells are treated carefully. That 
method also allows for a blocking of non-compact gauge fields, which is 



a{p) = d{—p)a{j))d{p) 



Another method, where different tastes contribute to a block variable, has been 
applied recently [54] to study the fourth root approach (cf. Subsection 3.4). That 
RGT drives the rooted staggered fermion to a sensible perfect action. However, the 
same is true for instance for the SLAG fermion [55], although the latter is incorrect 
under gauge interaction |i25j . as we mentioned before (in Subsection 3.4). 
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consistent in the sense that the hnk variables never connect fermionic 
variables on the same sublattice [57j . 

4.3 Application to the Gross-Neveu model 

We return to the Gross-Neveu model that we described previously in 
Subsection 3.5. More precisely we now consider its lattice formulation 
in terms of staggered fermions. Again we replace the 4-Fermi term by 
a Yukawa couplinj^ to an auxiliary scalar field 0. Since is taste-free, 
it is adequate to put its lattice variables on the cell centres z of the 
fermionic lattice [SB]. The standard formulation then couples 0^ in the 
same manner to the 2*^ taste variables located on the corners of the cell 
with centre z. 

As in Subsection 3.5 we considered the large limit, where the field 
<\)z freezes to a constant 0*^°^ Then the fermions can be integrated out, 
so that the RGT can be computed explicitly. In Ref. [SB] we derived the 
perfect staggered fermion action for this case. To analyse the scaling be- 
haviour, we evaluated two quantities of dimension mass for the staggered 
standard action and for the perfect action: 

• First we computed the chiral condensate (xx)- For the perfect 
action S\x, 0] this was achieved by a perturbation 

S\X. X, 0] = X, 0] + X, 0] . (4.14) 

The operator X has the standard lattice form X^xXxXx, and its 
perfect lattice form was computed again by the RGT technique, 
i.e. this perturbation was included to 0(e) in the transformation. 

• From the gap equation (analogous to the continuum eq. (13.371) ) 

20(0) = —^ I c/V In detM(p,0(°)) (4.15) 
(27r)^ Jb 

we extracted the fermion mass m j , which is dynamically generated 
by the breaking of the discrete, remnant chiral symmetry. M is 
the fermion determinant (see eq. fl3.18p ). either for the standard 
formulation or for the perfect formulation. 

In this context, we also considered the asymptotic scaling by in- 
vestigating how closely (f)^^\l/g) follows an exponential behaviour. 
(This behaviour is known in the continuum version of this model, 

^^By a Yukawa term we mean a product of a bosonic field and fermionic fields ^p, 
ip that contributes to the Lagrangian, as it also appears in the Standard Model. 
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see eq. 03.381) . and it characterises asymptotic freedom.) Theo- 
retically, asymptotic scaling does not need to be improved by the 
perfect action, since it is in principle independent from the scaling 
itself. Nevertheless we observed that it is significantly improved 
as well [56], in agreement with similar observations for truncated 
classically perfect actions for SU (3) gauge theory 



While these calculations involve lengthy expressions, the outcome for 
the (dimensionless) ratio of these two terms, which represents our scaling 
quantity, takes a simple form. 



ixx) _ 2sinh(am//2) (xx) 



standard aUl f 171 f 



= 1 . (4.16) 

perfect 



Hence the perfect scaling is indeed confirmed, i.e. for the perfect action 
the considered scaling ratio takes the exact continuum value at any lattice 
spacing a. In contrast, for the standard action this ratio is only obtained 
in the limit a — 0. We add that in this case also the classically perfect 
action scales perfectly; artifacts are switched off by the large N limit [56] . 



4.4 Exact supersymmetry on the lattice 

Since the RGT technique enables us to transfer continuum properties to 
the lattice without any damage in the physical observables, this procedure 
can in principle also preserve exact supersymmetry (SUSY) on the lattice 
[BOj. This may appear surprising, because continuous SUSY seems to 
contradict the lattice structure. For a review which presents a variety 
of approaches to handle SUSY on the lattice we refer to Ref. [61], and 
examples for further efforts to construct exact lattice SUSY are collected 
in Refs. [62]. 

For an illustration of the perfect action treatment of SUSY, we con- 
sidered the simplest supersymmetric model [63j: it is given in o? = 2 by 
the Lagrangian 

C[iIj, (fl = ^/'tA^ + d^ifd^ip , (4.17) 

with a Major ana spinor ip and a neutral scalar field (f. The action is 
invariant under simultaneous transformations with 



dip = --f^d^ipe , 6ip = eip, (4.18) 

where £ is a two-component Grassmann vector. The SUSY transforma- 
tion generators form a closed algebra with the translation operators, 

[5i, S2](p = {ei-ff,e2 - e2-f^ei)di,ip . (4.19) 
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By blocking from the continuum we transfer this model to the unit lattice 
and arrive at 

- (4.20) 

where G is the perfect fermion propagator (14. 3p . and Gg is a perfect scalar 
propagator (an obvious generalisation of the form given in eq. fl2.1ip ). If 
we perform the SUSY transformations f l4.18p in the continuum, they are 
carried over to the blocked lattice fields, 

5^^ = -7^ / d^ip{u)du e , 5^^ = s I ^{u)du , (4.21) 

which — under these transformations — obey exact SUSY too. 

In particular, we may treat the term = 7^<^ as a continuum current. 
As a general prescription we block a continuum current to the lattice by 
integrating its flux over the face f^^^^ between adjacent lattice cells [46], 

J^.,. = I d^-'yj,{y) ■ (4.22) 

This blocking scheme is illustrated in Figure [8] on the right-hand side. 
The lattice divergence of the blocked current is then equal to the conti- 
nuum divergence integrated over the corresponding lattice cells, 

5Jx = X](^M,^+A/2 - J,i,x^fi/2) = / d'^yd^j^iy) . (4.23) 

In this way, the transformations ( ]4.2ip can be expressed solely in terms 
of lattice quantities, i.e. the lattice current and the blocked lattice field 
[60] . Accordingly, the algebraic relation (14. 19p is now precisely reflected 
in terms of lattice quantities as 

82]^ = {eiV^J^e2 - £2V^ J^£i) , (4.24) 

where is the standard (symmetric) lattice derivative. 

These properties also extend to the free 2d Wess-Zumino model, which 
involves an additional scalar field that balances the fermionic and bosonic 
degrees of freedom, and to the free 4d Wess-Zumino model, where further 
field components are added. 

In terms of classically perfect fields, the continuum SUSY transfor- 
mations can be carried over to the lattice also in the interacting case. 
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However, the explicit construction of the corresponding lattice terms is 
a challenging numerical project, which has not been carried out so far. 

Still the lattice is hostile by its nature towards SUSY. In order to 
simulate SUSY models nevertheless, also the discrete formulation on a 
fuzzy sphere — that we already mentioned in Subsection 3.5 — should be 
considered [01]. For corresponding simulations we refer to Refs. 
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Figure 8: A cartoon of the schemes that we used to block various quanti- 
ties from the continuum: matter fields are blocked to the lattice by inte- 
grating the continuum field in a lattice cell, with the convolution function 
n. The perfect propagator G is obtained by integrating all continuum 
propagators between points in the corresponding lattice cells, as shown on 
the left (the formulae are given in Subsections 2.2, ^.l and 4.2). In the 
centre we illustrate the blocking for gauge fields, to be discussed in de- 
tail in Subsection 5.3. Here we integrate all straight parallel transporters 
between continuum points, which have the same relative position in adja- 
cent lattice cells. At last, a perfect current is obtained by integrating the 
continuum flux through the face between adjacent lattice cells ( it will be 
used again Subsection 5.3). 



5 Perfect Lattice Perturbation Theory 

On the level of analytical calculations, the construction of perfect lattice 
actions can be extended from the free fields to perturbative interactions. 
As a first example, we already sketched the computation of a perfect 
chiral condensate in Subsection 4.3. 

The method of blocking from the continuum is still applicable and 
highly efficient for this purpose. One now blocks various fields in such a 
way that all the continuum propagators between the continuum points 
in the lattice cells are integrated over. In the case of gauge interactions, 
also the gauge fields undergo a blocking procedure, which can be made 
explicit most conveniently for non-compact gauge fields, see Subsections 
5.3 and 5.4, and for illustrations Figures [8] and [TOl 
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5.1 The anharmonic oscillator 



As a toy model from quantum mechanics, we considered the anharmonic 
oscillator [67|. We write its action in field theoretic notation as 



s[^] = jdt[ + + \^{tY ] , (^(t) G IR) . (5.1) 

As in the case of the quantum rotor (discussed in Refs. [HI |18] and 
reviewed in Subsection 2.3), we use the ratio between the first two en- 
ergy gaps, Ai?i = El — Eq and Ai?2 = E2 — Eq, a.s a. scaling quantity. 
In continuum perturbation theory, the corresponding expansion can be 
found at many places in the literature, e.g. in Ref. [HH]- In terms of the 
dimensionless interaction parameter A := \/m? one obtains 

||(,) . 2 ^ 3A - if . !flA3 - 15|25£x- , O(A^) . ,5.2) 

First we evaluated this ratio to a high precision by Metropolis Monte 
Carlo simulations and we compare it to the perturbative results in various 
orders in Figure EI] (above). The latter approach the correct result only 
laboriously in a small range for A, even if we include the fourth order. 
This may serve as a caution to be careful in general with extrapolations 
to finite interaction strength based on perturbation theory. 

Next we calculated the perfect lattice action to 0(A). We chose the 
RGT parameter so that the action at A = consists of nearest neigh- 
bour couplings only. This is possible in Id field theory (i.e. quantum 
mechanics) with the parameter given in eq. (12.131) . We then extended 
the blocking from the continuum to 0(A). This generates additional 2- 
spin and 4-spin terms, which were written down explicitly in momentum 
space [HZ] . There inverse Fourier transform yields a set of couplings that 
we computed for various parameters rn? up to a coupling distance of 
two lattice spacings. This truncation is justified because the couplings 
undergo a fast decay, which speeds up for increasing w?. 

Finally we simulated the resulting perturbatively perfect action. For 
the scaling test, we fixed A = 0.005, i.e. a value where Figure ED (above) 
suggests the validity of first order perturbation theory. The results at 
various correlation lengths are compared to the outcome with the stan- 
dard action in Figure [5TT] (below). At a correlation length of ^ = 5 (in 
lattice units), both actions perform very well, and below 2.5 both suf- 
fer from similar scaling artifacts. In between, there appears a window 
where the perturbatively perfect action seems superior, as we observed 
at ^ = 10/3. 
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Figure 9: Above: The ratio between the leading energy gaps, AE2/AE1, 
against the perturbative predictions. Higher orders do extend the range 
of the validity of perturbation theory gradually, but even at fourth order 
it is still limited to small interaction parameters (X <0.03). 
In the plot below we compare simulation results at X = 0.005 with the 
standard action and the 0(A) perfect action, for different correlation 
lengths ^ in lattice units fUj 



5.2 The Yukawa term 

We computed a perturbatively perfect action in the framework of the 
Gross-Neveu model with staggered fermions, cf. Subsections 4.2 and 4.3, 
but now for four tastes. In this case (without a large N limit), the 
auxiliary scalar field $ is not constant anymore, but we assumed it to 
be small. More precisely, we absorbed the Yukawa coupling in $ and 
considered its first order. In this approximation — which describes the 
asymptotically free system at high energy — the perfect staggered action 
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takes the form 

xy,ij z 

+ Yl -z,y- ^)xi^z , (5.3) 

xyz,ij 

where x, y run over the lattice which hosts the fermionic degrees of free- 
dom, whereas z runs over the plaquette centres, and i,j = 1...4. If we 
take the spacing between fermion components of the same taste as the 
unit, z is spaced by 1/2. In momentum space we write the interaction 
term as x{—p)o'{p, q)x{~(l)^{p + Q)- the taste space, the shifted kernel 
5" is a 4 X 4 matrix, which only couples tastes of the same sublattice. To 
be explicit, its first order perturbatively perfect form reads [SB] (we use 
the notation of eq. (I4.13P ) 



a{p,q) = d{p)a{p)a{p,q)a{-q)d{-q) 



with the matrix elements 

{l,n) (m,n) 











-0"3 \ 







^3 








-CTg 







\ 5-3 








-^0 / 



(5.4) 



rniljTi) 2 «(m,n) 2 

^i^+m,+n^ PuqAp^q)u 

(/,n) (m,n)r (Z,n) 

u=i Pf qi^ [P'^ ^ qt^ 



{l,n) (m,n) 

C3{P,q) — p(«,n)2^(m,n)2 



X 



11'^ (/,n) (m,n)r (l,n) , (m,n)i ' ^'-^''^^ 

u=i Py qv \Pv + qu \ 



where p^]" = p^ + 47r/^ + 27m^, and G {0, 1}. Note that these matrix 
elements are 47r periodic, in accordance with the central positions of the 
auxiliary scalar variables. 

This Yukawa term identifies the direction of a "renormalised trajec- 
tory" (a line ofperfect actions in parameter space) emanating from the 
critical surfaceo The corresponding couplings in coordinate space can 
be evaluated numerically, and they have been applied — in a truncated 
form — in lattice simulations 1691. 



^^The endpoint of this trajectory is the free perfect action that we identified before, 
due to the asymptotic freedom of the Gross-Neveu model. 
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5.3 Perfect gauge actions and the axial anomaly 

The attempts to formulate non-local lattice fermions with a finite gap 
at the edge of the Brillouin zone were unsuccessful; we mentioned the 
SLAC fermion in Subsection 3.4. A refined approach was presented by 
C. Rebbi, who formulated a non-local fermion with divergences at these 
edges instead ^70j. However, the Rebbi fermion does not reproduce a 
non-zero axial anomaly, as A. Pelissetto pointed out [71j. 

If we construct the perfect fermion for a 5-function blocking RGT, 
we obtain a non-locality of the same type as the Rebbi fermion [42l H3] , 
hence we wondered what happens to the axial anomaly in that case. 

Once we deal with the (5-function blocking, the relations fl4.ip turn 
into equations. In momentum space they read 

i|r(p) = ^ n(j9 + 27r/)^(p) , ^{p) = ^ n(p + 27rl)ip{p) . (5.6) 

Analogously, we now block an Abelian gauge field Afj,{x) from the con- 
tinuum to construct the non-compact link variable [461 fT5] 



A^,x= j d^u{l + u^-x^)A^{u)+ i (fv{l-v^ + x^)A^{v) . 

(5.7) 

Here x is a link centre on a unit lattice, so that we integrate over adjacent 
lattice cells. This blocking scheme is illustrated in the centre of Figure 
[HI A gauge transformation in the continuum, Ap, — >■ Ap, + 9^A, induces 
exactly a lattice gauge transformation 

Ap^^ Ap^^ + A^+^/2 - A^_^/2 , A^ = / (X^y \{y) , (5.8) 



which shows that this lattice gauge field is covariant. In momentum space 
it takes the form 

Ap{v) = J2 n^lP + 2nl)i-iy^ A,{p + 2nl) , 

where H^(p):=^H(p) (5.9) 

Ptj. 

is anti-periodic over the Brillouin zone in the /x-direction (and periodic in 
the other directions). It is convenient to start from a continuum action 
in the Landau gauge, which leads to the perfect lattice gauge action |15] 

S[A] = -^^j^d''p\A,{-p)Gf{p)-'A,{p), 

lp{p + 2nl 
{p + 2Txlf 
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where the term i?^^-* smears the RGT in analogy to the treatment of 
the matter fields in eqs. fl2.5l) and fl4.2l) . As a remarkable property, the 
specific choice 

R^'\p) = 1-^pI (5.11) 

yields for an Abelian gauge field in d = 2 the standard plaquette action, 
which is therefore perfect already [T5]. This property is similar to the 
standard lattice scalar and the Wilson fermion in d = 1, which are also 
perfect in the non- interacting the choice of suitable RGT para- 

meters confirms (this was pointed out previously in Subsections 2.2 and 
4.1). Again we use this property as a tool to optimise the RGT in view 
of locality also in higher dimensions. 

In Ref. |16] we considered a more general perfect lattice gauge action, 
without previous gauge fixing in the continuum. There we also proceeded 
to perturbation theory to the first order in the gauge coupling g. In this 
case, the blocking RGT of the fermion field is extended from eq. (15. 6p to 
the form 

^,{p) = 5^ n(p + 27rl)Up + 2vr/) + [ 

X K^{p + 27r/, q + 27r0^^(p - q)Xij^P,{q + 27rl) , (5.12) 

and correspondingly for ^'j. Eq. (15.121) refers to a SU{N) (or U{N)) 
gauge field, where A*^ are Hermitian generators. However, in the current 
Subsection we will deal with the Abelian gauge field, where the last factor 
simplifies to Af^{p — q)ip{q + 2txI). 

The kernel has to be regular, and gauge covariance requires 

(Pm - q,)K,{p, q) = Tl{p - g)n(g) - U{p) . (5.13) 

To this order, the perfect action includes — in addition to the pure 
fermion and pure gauge action — an interaction term of the structure 

V[^>,m,A] = l^^d'^pd''q^,i-p)V,ip,q)Alip- q)X^^^,iq) . 

(5.14) 

The explicit form of Vfj_{p, q) is rather lengthy; it is written down for the 
gauge group U{1) in Ref. [16], and for QCD in Ref. [153. Gauge invariance 
requires that the vertex function obeys the lattice Ward identity 

if^q)Mp, q) = Giq)-' - G{p)-^ , (5.15) 

where G is the perfect fermion propagator (14.31) . 
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To investigate the axial anomaly, we now block the axial continuum 
current 

to the lattice. Following the prescription fl4.22p we obtain a lattice current 
J^j,. In momentum space it takes the form 

In the Schwinger model (two dimensional QED [72]), the continuum cur- 
rent in a gauge background is known to obey the relation 

{jlip))A=-^e,pP.A,{p) , (5.18) 

which can be derived for instance with dimensional regularisation. We 
integrate out the continuum gauge field and take the lattice divergence 
of the current J^^ in the perfect lattice background, which leads to 

i^Mp))-* = f E'-7r^|^n(P + 2.;)(-i)'' 

X n,(p + 27r/)G^^(p)-iA.(p) , (5.19) 

where G^p^{p) is the perfect gauge propagator (without previous gauge 
fixing in the continuum it generalises to a tensor). On the other hand, 
we build the perfect topological charge density in agreement with the 
blocking recipe of Subsection 2.3, i.e. we block the continuum density 
^^^J.u^|J.Al, to the lattice. We find that it coincides precisely with the 
lattice divergence of the perfect current [16], 
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:=- d'ye^^d^AM = {^JDa ■ (5.20) 

J Cx 

Therefore the perfect lattice action does indeed reproduce the perturba- 
tive axial anomaly correctly, at any lattice spacing. 

Hence the perfect fermion constructed with a (5-function RGT is the 
only non-local lattice fermion that did not run into conceptual trouble. 
It represents the only conceptually successful implementation of chiral 
symmetry in its standard form on the lattice. For practical purposes, 
however, one prefers the local form, which is still compatible with a 
modified chiral symmetry, as we will discuss in Section 7. 
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5.4 The perfect quark gluon vertex function 



The kernel function in the fermionic blocking scheme to the first order 
in the gauge coupling is submit to the constraint fl5.13p . and it is difficult 
to find explicit solutions for it. This issue is discussed in Ref. [50]. In the 
small momentum expansion one obtains unambiguously to the leading 
orders 

K^{p,q) ^ ll [i + Y^{p'-4[P;.(PM-g;.)+gJ] -5[p{p-q)+f]^ 

(5.21) 

where p = {p^,p), q = {q^,q)- We also discovered two full solutions, 
which can be evaluated numerically. We display here one of them, de- 
noted as the "recursive kernel", in a corrected form (unfortunately this 
formula contains an error in Ref. [50]). 

^.(^'^) = iEn|^n,((p-g)/2")n(g/2-^^) 
X sin(g^/2"+i)/C^(p/2"+\g/2"+^) 

MP,?) = E [n'"'^^- + "^'-^'^l"f-""^'-^'^] , (5.22) 
re{o,i}''-i ^^^^ (1 + ^ ) 

where / excludes the /i-component. 

Along with eq. (15.14^ this provides a fully explicit — though some- 
what complicated — form of the perfect quark gluon vertex function [TB] . 
We recall that its gauge invariance is guaranteed by the Ward identity 

We worked out a truncated version of the quark gluon vertex function 
in coordinate space, which we applied in simulations to be addressed 
in Subsection 6.4. Its general form involves terms in the full Clifford 
algebra. It simplifies drastically if we map the system down to c? = 2 by 
dimensional reduction. For that case we gave explicit couplings, including 
the perturbatively perfect clover term, at various fermion masses in Ref. 

m. 



6 The Hypercube Fermion 

6.1 Construction of the Hypercube Fermion 

In Subsection 4.1 we have already described the truncation of the perfect 
free fermion to a hypercube fermion (HF) by means of periodic boundary 
conditions. We gave the couplings for the massless fermion in Table [TJ 
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A^{p - q) 




Figure 10: Overview over the construction scheme of the perfect quark 
gluon vertex function: the lattice fermion fields ^ and \1/ and the gluon 
field are obtained by blocking from the continuum, cf. Figurel^ For a 
given point (as the one marked by the central dot in this Figure) we then 
integrate the continuum propagators to all continuum fields involved in 
this blocking process. 

We have also seen that it has excellent scaling properties, and we will 
illustrate in Section 7 that its approximation to chirality is excellent as 
well. 

The numerical treatment of all couplings inside a unit hypercube is 
clearly more complicated than the effort for standard formulations (Wil- 
son or staggered), which only deal with nearest neighbour couplings. 
However, simulations with this generalised form are feasible and have 
been carried out. To this end, the first question was how to handle the 
link variables to include the gauge interaction. If the spinors and 
are coupled in some lattice action, we can arrange for gauge invariance of 
the corresponding term by connecting these spinors over a lattice path, 
where the compact link variables are multiplied [28j (cf. Subsection 3.4). 
The way to do so is ambiguous. The perfect or classically perfect actions 
do actually determine the couplings to these lattice paths (once the RGT 
is chosen). But their evaluation is tedious, and a truncation of the path 
length is required, which is again arbitrary. 

The simplest approach just includes connections over the shortest 
lattice paths. For most connections in a hypercube there are several 
shortest paths, which are then all included with the same weight and 
averaged over. In this way, we introduce "hyperlinks" which connect a 
lattice site x to all the 3^^ sites in a common unit hypercube with x. 

(2) (3) 

An illustration of 2d and 3d hyperlinks, f/^+i, and U^^^j^p , is given in 
Figure [TTl Regarding the numerical implementation, it is favourable to 
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Figure 11: A 3d illustration for the construction of hyperlinks in the 
gauging of a hypercube fermion. 



construct these hyperlinks in a given gauge configuration recursively, i.e. 
to start with the 2d hyperlinks (over plaquette diagonals), which then 
also serve as building blocks for the 3d hyperlinks, from which one arrives 
conveniently at the 4d hyperlinks [73j . 

In Ref. [73] we discussed a suitable preconditioning method for the 
HF. The goal is to divide the lattice into sublattices such that a block 
structure in the lattice Dirac operator Dx,y{U) is generated. This block 
structure allows for the transition to an equivalent Dirac operator with a 
better condition number. This means that the ratio between the maximal 
and minimal eigenvalue of D'^D decreases, and this ratio is essential for 
the computational effort in a simulation. Moreover, the transformed 
operator is block diagonal, which also simplifies its inversion. 

For the Wilson fermion, this method is usually apphed with two sub- 
lattices, which do not contain nearest neighbour sites, as distinguished 
by the sign term fl3.30p . It is known as even-odd (or red-black) precondi- 
tioning. For the HF a set of 2'^ sublattices is suitable, which we denoted 
as "rainbow preconditioning"]^ In fact, it yields gain factors between 
3 and 4 in typical QCD simulations [73]. (We add that also for Wilson 
fermions it can be profitable to deal with a larger set of sublattices, as 
the locally-lexicographic symmetric successive overrelaxed preconditioner 
(//-SSOR) shows [B].) 

6.2 Approximate rotation symmetry 

A first simulation of this HF with the simple gauging described above 
was presented in Ref. [SO]- We set the bare fermion mass to zero and 

^^This approach was also discussed for the hypercube scalar in Ref. [6]. 



54 




Figure 12: The "speed of light" in eq. h6. 1\) determined from the disper- 
sion relation of a pseudoscalar meson on the lattice. We see that the 
continuum relation Cmeson = ^ is approximated very well for the HF, but 
not for the Wilson fermion [5U^ . The simulation was performed quenched 
with the Wilson gauge action at (3 = 5. 



evaluated the dispersion relation for the pseudoscalar meson with the 
Wilson plaquette gauge action (see e.g. Refs. [2S]), at /5 = 5 in a quenched 
simulation on a lattice of size 6^ x 18. The corresponding "speed of light" 



_ M2 
\P\ 



(6.1) 



is shown in Figure [12] and compared to the result for the Wilson fermion. 
In this formula, E, M and p are the energy, mass and 3-momentum of 
the pseudoscalar meson. The continuum behaviour, Cmcson = 1, is marked 
by a dotted line. We see that it is approximated very well for the HF 
— which leads to Cmeson = 1.04(5) — in striking contrast to the Wil- 
son fermion. This property corresponds to an excellent approximation 
to Lorentz symmetry (resp. to Euclidean rotation invariance). However, 
the meson mass is strongly renormalised in this case. In lattice units it 
amounts to M ~ 3, hence it can hardly describe a pion. Of course, we 
were dealing with a very coarse lattice. Still, this property calls for a 
closer look at chiral symmetry, which we will undertake in Section 7. 



The HF has been applied successfully in simulations of the Schwinger 
model [72] with two flavours, on a 16 x 16 lattice at /3 = 6, 4 and 2 [75] . 
In these simulations, the gauge configurations where generated quenched. 



55 



_ 1e-04 
1e-05 
1e-06 



t 




o Wilson 
• fixed point 
OTP-HF 


o 




. SO-HF 








- ° o 












fib 













|x| 



Figure 13: The isotropic correlation function 1^6. S\} for the Wilson 
fermion, a fixed point fermion and two HF versions in the Schwinger 
model. We see in all cases but the Wilson fermion a smooth decay, which 
confirms the good quality of approximate rotation symmetry down to short 
distances. 

but the measurements did include the fermion determinant @ We used 
the Wilson plaquette gauge action, which is perfect for pure 2d U{1) 
gauge theory (see Subsection 5.3). 

First we present another test of the quality of rotation symmetry. 
Figure [13] shows the decay of the correlation function 

C3(x) = (foc^s^o-'S'.as^x) (6.2) 

against the distance where x are lattice sites in all directions. We 
measured this correlator for the Wilson fermion and for several types of 2d 
HFs, which are described in Ref. [7S]. For the TP-HF, the fermionic coup- 
lings correspond exactly to the description in Subsection 4.1, whereas the 
SO-HF is further optimised with respect to the scaling behaviour. Both 
variants also include a clover term. This plot shows in addition the results 
obtained with the classically perfect action constructed in Ref. |13j . 

The observation that the HFs and the (far more complicated) classi- 
cally perfect action display a much smoother decay of C3 down to short 
distances |x| approves their good approximation to rotation symmetry. 

^°In the recent literature, it became fashionable to denote this kind of simulation 
simply as "dynamical" , although the fermion determinant is still treated as constant 
in the generation of the configurations. 
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Figure 14: The dispersion relations for the "pion" and the "r]-meson" 
in the Schwinger model with various lattice fermion formulations: Wil- 
son fermions (diamonds), the classically perfect action (filled circles) and 
three types of HFs, in particular the scaling optimised SO-HF (little empty 
boxes) /yj/ . The solid lines mark the behaviour in the continuum. 

6.3 Spectral properties in the two-flavour Schwinger 
model 

6.3.1 Applications of the hypercube fermion 

In Ref. [75] we also measured the dispersion relations for the states, which 
are analogous to the pion and the 77-meson (a general discussion of these 
properties of the Schwinger model can be found in Ref. [76]). Again, for 
increasing momenta (in lattice units) scaling errors due to lattice artifacts 
are enhanced, cf Subsections 2.2 and 4.1. In Figure we show these 
dispersions, which are again strongly improved for the classically perfect 
action and for the HFs, in particular for the SO-HF. It is remarkable 
that the latter — which is still very simple — scales at least as well as 
the highly involved classically perfect action of Ref. [13] . 

6.3.2 Applications of a truncated perfect staggered fermion 

In Ref |77| we constructed with a similar procedure a truncated perfect 
staggered fermion, starting from the formulation described in Subsection 
4.2. We applied it in Schwinger model simulations as well, and these 
simulations were truly dynamical. We designed a variant of the Hybrid 
Monte Carlo algorithm [78] , which is particularly suitable for this formu- 
lation. It uses a simplified action (the standard staggered fermion action 
with fat links) for the Molecular Dynamics, and the full quasi-perfect 
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action in the Metropolis accept/reject step. This provided a numerically 
cheap evaluation of the force, along with a still useful acceptance rate 
(as expected, the latter decreases at strong gauge coupling, which corre- 
sponds to a large physical volume). Again we used a 16 x 16 lattice and 
the plaquette gauge action. At /3 = 3 we found the neat vr and rj dis- 
persion relations shown in Figure [T^ In the framework of that project. 
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Figure 15: The dispersion relations for the "pion" and the "rj-meson" in 
the 2-flavour Schwinger model with dynamical staggered fermions. We 
show results for standard staggered fermions and for truncated perfect 
staggered fermions /77/, similar to the HF discussed before. The solid 
lines mark the behaviour in the continuum, which is also here very close 
to the data for the truncated perfect action. 



we further investigated the "meson" masses under variation of the gauge 
coupling (3 (the parameters are given in detail in Ref. [T7j). The results, 
shown in Figure [161 reveal again that the truncated perfect formulation 
is far more suitable to approximate the continuum behaviour, and in 
particular to realise light pions — even on coarse lattices. 

This work has pioneered on one side the construction and application 
of improved staggered fermions — which became indeed fashion in the 
beginning of this century — as well as the use of a simplified force term in 
the Hybrid Monte Carlo simulation of a complicated quasi-perfect action. 

6.4 The charmonium spectrum 

Regarding QCD, we performed simulations to evaluate the charmonium 
spectrum employing the HF [79]. In this case we used a truncated version 
of the perfect quark gluon vertex function discussed in Subsection 5.4. 
The result is shown in Figure [T71 This was a quenched simulation with 
the Wilson gauge action at /3 = 5.6 on a 8^ x 16 lattice. The bare quark 
mass was set to m = 0.9 — that value was adapted to match the ground 
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Figure 16: The "meson" masses in the Schwinger model from dynamical 
simulations with different types of staggered fermions, at various values 
of P, resp. lattice spacings a oc We confirm that the results for 

the truncated perfect staggered fermion are far closer to the continuum 
results, and they provide much lighter "pions". 



state ?7c(2.98GeV). Considering the fact that only this ground state was 
used as an input, the agreement with the experimental spectrum is clearly 
satisfactory. The inclusion of a term oc •y^'julp in the vertex function (see 
also Ref. [50]) was especially helpful for the quality of the charmonium 
spectrum. 

GeV -- 



3.5 — 



Vc J/i5j 

Figure 17: The charmonium spectrum, measured in simulations with the 
HF and a truncated perfect quark gluon vertex function l79^ . The exper- 
imental values are dashed, and the ground state of t]c sets the scale. 



6.5 Spectral functions at finite temperature 

At last, we add that the HF is currently being applied in studies of the 
spectral functions of QCD at finite temperature pniEl]- These spectral 
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functions, depending on the frequencies, are obtained from lattice data 
using the Maximum Entropy Method, which was suggested for this pur- 
pose in Ref . [82] . For the HF they reveal a continuum-like behaviour up 
to much larger energies than it is the case for the Wilson fermion, see 
Figure [m The basis for this virtue is that the HF moves all the doublers 
to a unique cutoff scale, whereas the Wilson fermion splits them into 
four (in general d) subsets, see Subsection 7.5. We recall that the naive 
doubler species have between 1 and d momentum components vr/a in the 
free case. This causes the splitting in their cutoff energy for — an 
effect that Drf overcomes. 
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Figure 18: The spectral function aps, as a function of the frequencies uj, 
at critical temperature Tc = oo (on the left, free case) and at finite Tc 
(on the right, interacting case). These results are obtained from lattice 
data using the Maximum Entropy Method for the HF and for the Wilson 
fermion. 



7 Chiral Correction by means of the Over- 
lap Formula 

7.1 The Ginsparg- Wilson relation 

In Subsection 4.1 we discussed the block variable RGT for the free fermion 
field. Eq. (14. 2 p introduced the transformation term, which we assume (for 
simplicity) to be diagonal in the spinor indices. 
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For the bare fermion mass m = the inverse perfect lattice Dirac oper- 
ator has the structure 



Rxy — „^xy i^-^) 



+ ,7,2) 

where we are using the notation of eq. (14 .7^ and = p^p^ (in coordinate 
space, the products are lattice convolutions). We repeat that in the limit 
of a 5-function RGT, — > 0, the lattice action is invariant under the stan- 
dard chiral transformation fl3.14p . as we see from the anti-commutator 
{^,75} = 0. 

However, we also saw that locality requires 7^ 0, in agreement with 
the Nielsen Ninomiya No-Go Theorem [22], and we found for m = the 
term 

1 

— ( 
2 

to be optimal for locality [15]. In this case, we have to face {D,75} 7^ 0, 
but due to our requirements for R in Subsection 4.1 (eq. (14. 5p ) the anti- 
commutator 

{^-1,75} = 275/2 (7.4) 

is local. The exact form of the factor 11 (which depends on the blocking 
scheme, for our case it is given in eqs. fl4.3p and fl2.1ip ) does not affect 
this relation. This already indicates that for a given term R there is a 
variety of solutions to eq. (17. 4p . 

We stressed in Section 4 — and a long time ago, for instance in Ref. 
[56] — that the specific chiral symmetry breaking by the term R cannot 
distort any physical properties, in particular not those related to chirality. 
Relation (17.41) illustrates this again, since a local term R cannot shift the 
poles in the propagator D^^ . If we multiply the operator D from both 
sides, we arrive at the equivalent equation 

{D,-i^} = 2D-i^RD , (7.5) 

which is now known as the Ginsparg-Wilson relation (GWR)@ The di- 
mensions reveal that the right-hand side is 0(a) suppressed. As we men- 
tioned earlier, the GWR was written down in Ref. [44j, which was for- 
gotten over a long period until its re-discovery in 1997 |45j . 

We may compare this property to the Wilson operator, which breaks 
chiral symmetry such that {-Dw,m=o,75} amounts to a local term, unlike 
relation (17.40 . If we try to insert the free operator i5w,m=o into the GWR, 
we arrive at a non-local term R^y, which decays only as |x — ?/|~^ in d = 4, 



■^^If we allow a general Clifford structure for the term R, the right-hand-side turns 
into D{-f5,R}D. 
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and as \x — y\ ^ in d = 2 [H3]- Hence the Wilson operator does not solve 
the GWR. 

We continue to use the GWR as a requirement also in the interacting 
case. For gauge interactions, D{U) involves the link variables, and so 
does i?, if it couples different lattice sites. However, we stay now with 
the choice (17.31) . for which the GWR takes the simple form 



R.V = ^p6^y,p = l {D,75} = D^,D . (7.6) 



M. Liischer observed that this relation assures the exact invariance 
under a lattice modified chiral symmetry transformation [IH]. In the 
notation of eq. (13.141) . it is sufficient to consider this modified transfor- 
mation to 0{a) (due to the Lie group structure of the chiral rotation). 



\1/Z)\1/ + a\E' 



{D,j,}-Dj,D ^ + 0{a^) . (7.7) 



The GWR corresponds exactly to a vanishing term in the square bracket. 
In fact, it could have been also discovered by requiring invariance under 
such a lattice modified chiral symmetryjfl Unlike the remnant chiral sym- 
metry of staggered fermions, this chiral rotation involves the full number 
of generators that appear in the continuum theory. The continuum limit 
a — >■ yields a smooth transition to the standard form of chiral invari- 
ance. Based on these observations, even a non-perturbative formulation 
of chiral U{1) gauge theory (cf. Subsections 3.3 and 3.4) has been worked 
out on the conceptual level [HI]. Earlier works in this direction, based on 
the overlap formalism, are collected in Ref. [85] . 

Perfect Dirac operators solve the GWR, but — as we saw — they can 
in general not be constructed explicitly. A step towards an applicable 
solution was the observation that also classically perfect Dirac operators 
(cf. Subsection 2.3) are Ginsparg- Wilson operators [86] . Still, a trunca- 
tion of the couplings and therefore a deviation from exact chiral symme- 
try is needed in its construction, but in view of the possibilities to build 
approximate solutions this is a more accessible starting point than the 
(quantum) perfect action. 



^^The GWR is presented along these lines in Ref. [83], which also keeps a general 
term R in the transformation (17. 7p . Then the transformation of D reads 
D^[l + a{l- DR)-f5]D[l + 075(1 - RD)] =D + a{{D, 75} - D{R, ro}D) + Oia^). 
Again the invariance to 0{a) is equivalent to the GWR, and for local terms R there 
is a smooth transition to the standard chiral symmetry in the continuum limit. 
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A fully explicit solution was presented by H. Neuberger [SZ]. Let us 
recall the properties (13.161) and start from some massless lattice Dirac 
operator Dq, which we assume to be 75-Hermitian, i.e. 

Dl = 75^075 . (7.8) 

This property holds for practically all Dirac operators that have been 
considered, in particular for the Wilson operator -Dw@ Now we define 
the shifted operator 

A:=Do-l, (7.9) 
which is unitary if Dq is a GW operator. 
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Do75^o-{^o,75} + 75j • (7.10) 

Of course, we cannot assume this for a quite general Dq, and we pointed 
out before that it does not hold in particular for But we can trans- 
form A such that unitarity is enforcedly 

A Aov := A/VaU AlA,, = 1 . (7.11) 

In this way, we obtain the overlap Dirac operator 



L>ov = l + A, 



l + iDo-l)/^JiDl-l)iDQ 



= 1 + 75^ , H:=^,A, (7.12) 

where H is Hermitian, H = W . H. Neuberger introduced this operator 
[87] with Dq = D^^, and we denote the corresponding overlap operator 
as the Neuberger operator D^. 

Since the resulting lattice action S[^, ^E', U] has now a modified but 
exact chiral symmetry, one may be worried about the fate of the axial 
anomaly. However, the anomaly is in fact captured correctly, due to the 
variance of the functional measure under modified chiral rotations ^9] , 
which is analogous to the continuum. This property means an explicit 
realisation of the program formulated in Ref. |89]. The axial anomaly 
has been computed perturbatively by many authors, for instance in Refs. 
[m HHl EEl [90] for general Ginsparg- Wilson operators, and specifically 
for also in Refs. [91]. A proof which extends to all topological sectors 
was given for in Ref. [92] • In fact, this extension is non-trivial, as the 
considerations in Refs. [93l underline. 



An exception is the operator of the so-called "twisted mass fermion" iSS] . For the 
staggered fermion, the analogous relation e{x)pfj_{x,y)e{y) = — /9^(x,?/) holds, where 
e{x) is the sign factor (|3.30[) . and is the nearest- neighbour coupling in eq. (|3.29p . 

^''Note that A does in general not commute with A^A, so we have to specify an 
order where to multiply the inverse square root. This is not necessary anymore in the 
form given in the lower line of eq. ()7.12|) . 
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7.2 Massless lattice fermions in d = 3 



In three dimensions, there is no chiral symmetry (since there is no matrix 
75 at hand), but parity takes a role similar to a discrete chiral symmetry!^ 
The parity operator TZ is then the analogue to 75, and we can write the 
parity transformation on the lattice for spinors and compact link variables 
as 

P : i^^n , ^ ^7^^, , 

U,,. - U^,^ := Ul_^ . (7.13) 

Similar to the 75-Hermiticity (17. 8p . practically all lattice Dirac operators 
which are considered obey 

D{Uy = nD{U^)n . (7.14) 

The relation 

D{U)+D{U)^ = Q (7.15) 

then implies parity invariance of the action 5* = 'Ylixy^ ^^i^)xy^ ^^"^ 
fermion mass zero. However, this condition (17.151) is again not easy to 
fulfil — it leads to a doubling problem as the chiral symmetry does in 
even dimensions. The doublers in the naive action can be avoided by 
a Wilson term, but this term breaks the condition (17.151) and therefore 
parity symmetry, so that additive mass renormalisation sets in. 

A (massless) overlap fermion in ci = 3 was introduced in Ref. [91] , 
which was actually the first place where the overlap formula (I7.12p ap- 
peared. In Ref. [95] we considered generally a 3d Ginsparg- Wilson oper- 
ator given by the condition 

D + D^ = D^D, (7.16) 

and we studied a lattice modified parity symmetry. This modification 
alters the transformations (17.131) such that 

^,^tn{l~D)^, . (7.17) 

For a solution D to eq. (17. 16p the lattice action is exactly invariant under 
this modified parity symmetry, but the functional measure transforms as 

V^V^ [<lei{l ~ D)]-^V^! Vm . (7.18) 



■^^We encountered a discrete chiral symmetry before in the Gross-Neveu model (in 
the continuum), see eq. p.34p where in view of the chiral rotation (|3.14p only the 
angles a mod 27r G {0,7r} occur. In that case, extending the 4-Fermi term to 
(^\I>)2 ^ j-^^^^^[;-j2 pj-Qijiotes the chiral symmetry to a continuous form (Thirring 
model). 
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Once more in analogy to the chiral symmetry in even dimensions, this 
transformation of the measure gives rise to the requested parity anomaly 
[95j . We should clarify that this is not an anomaly in the usual sense, 
which has a unique value, but it comes with an arbitrary integer factor (a 
comprehensive discussion is given in Ref. [96] )• Hence a successful regu- 
larisation should capture all possibilities for this integer. In the current 
case, they are all captured indeed by varying either the kernel R in the 
GWR (which then modifies the right-hand-side of eq. fl7.16p to 2D^RD), 
or by considering different kernels Dq in the overlap formula fl7.12p . 

At this point we add that also a pure Abelian 3d Chern-Simons gauge 
theory with the continuum action 

S[A] = j d'xe^.pA^.d.Ap (7.19) 

suffers from a doubling problem on the lattice due to the occurrence of 
a linear momentum [97]. A solution can be found also here either by the 
perfect action formalism, or by a formula of the overlap-type [98] . 

7.3 The overlap hypercube fermion 

We proceed to a slightly more general form of the term R in the GWR 
(17.51) : we now allow for a parameter p > 1 in eq. (17. 6p . The overlap 
formula can be adapted to general forms of R [991 111]- Fo^' the still 
simple form that we are considering now, it reads 

^ov = p(l+A/v^) =p(l+75^) , 

A ■= Do-p = ^5H. (7.20) 

As we mentioned, the standard formulation inserts as a kernel the Wilson 
operator, Dq = -Dw; which is far from chiral, and which undergoes a 
drastic change in the overlap formula to yield the Neuberger operator 
-Dn- Ref. [99] suggested to consider more general possibilities for Dq, and 
motivated in particular the use of a truncated perfect operator, which is 
approximately chiral already. In this case, the square root in eq. (17.20p 
is close to the constant p, and the transition from Dq to Dov is therefore 
only a modest modification. An intuitive argument for this property is 
that an exact GW kernel Dq reproduces itself identically in the overlap 
formula (I7.20p (for a fixed parameter p)@ 

^^In this sense, the overlap formula captures all solutions to the GWR. Attempts 
to construct chiral lattice Dirac operators beyond GWR solutions were considered in 
Refs. [Too] . 
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In particular we are using D^p as the input Dirac operator — note 
that it is 75-Hermitian as well. Its inexact chirality is then corrected by 
the overlap formula, which leads to -DqvHF, the operator of the overlap- 
HF, while keeping both Dirac operators similar, 

D0.RF ~ DuF . (7.21) 

The motivation for this construction is that the property (17.211) allows 
us to preserve other virtues of D^p (beyond the approximate chirality) 
in the chirally exact formulation DqvHF- As such virtues we are going to 
discuss: 

• a high level of locality 

• approximate rotation symmetry 

• a good scaling behaviour. 

Below we will summarise results for these three aspects and comment on 
their meaning^ Still before that we mention that the simulation of an 
overlap fermion with a hypercubic kernel requires more computational 
effort compared to -Dn- The use of the complicated kernel by itself cor- 
responds to an overhead of about a factor 15 in QCD (without applying 
the preconditioning method reviewed in Subsection 6.1). However we 
should consider that in 4d simulations, the inverse square root in the 
overlap formula fl7.20p (resp. the sign function -^^^ ) is approximated by 
polynomials. Again thanks to the relation (17.211) . the convergence in the 
polynomial evaluation is faster — for a fixed precision one gains back by 
this property about a factor of 3 lW2\ \103\ 1104] . Hence an overhead of 
about a factor 5 remains. We are confident that this factor will be more 
than compensated by the superior properties listed above, which we will 
discuss in Subsections 7.6 and 7.7. 



7.4 The axial anomaly in the continuum limit 

On the conceptual side, we were able to prove that also the overlap- 
HF (obtained by inserting Dq = Dhf in the overlap formula, as we just 
advocated) has the correct chiral anomaly in the continuum limit of all 
topological sectors |105j ,^ This proof required a number of generalisation 



^^Of course, one could also construct an approximate Ginsparg- Wilson operator 
directly by starting from some short-ranged parameterisation ansatz and tuning the 
couplings such that the GWR is minimally violated. This was done in Refs. ^101] for 
the Schwinger model and for QCD. 

''^This may be compared to the fully perfect action, which even displays the correct 
axial anomaly at finite lattice spacing, as we discussed in Subsection 5.3. 
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steps compared to the proof that had been worked out previously for the 
Neuberger operator [22] • 

In Ref. |105j we gave a non-perturbative evaluation of the continuum 
limit of the axial anomaly and the fermion index. To this end, we formu- 
lated the Dirac operator Dhf on a 2n-dimensional Euclidean lattice in 
a form, which is well-suited for analytic investigations. We used it first 
to study the dependence of the doubler structure of Dhf on its coupling 
parameters. Then we evaluated the classical continuum limits of the ax- 
ial anomaly and the index of the overlap-HF operator, showing that the 
correct continuum expressions are recovered when the parameters are in 
the physical (doubler-free) range. This continuum limit relies only on 
general properties of D^p and not on its explicit form (in contrast to the 
previous evaluations in the Wilson case |92]). The main new technical 
tool was a set of identities, which allowed the continuum form 



of the axial anomaly to be extracted, and its coefficient to be evaluated 
topologically. 

These properties are basically not specific to the HF structure. We 
expect that this proof can be extended to completely general overlap 
Dirac operators obtained by inserting an ultralocal (and 75-Hermitian) 
lattice Dirac operator Dq (involving the full Clifford algebra of 7 matrices, 
as the operators in Refs. [HI llOll I106j ) into the overlap formula. 

7.5 Approximate chirality of the hypercube fermion 

An exact solution to the GWR with a term R of the form (17. 5p and p > 1 
has its spectrum cr(D) on a circle in the complex plane with centre and 
radius p, as we see from the relation A'^A = p^. We call it the GW circle, 



This property reveals immediately the absence of additive mass renor- 
malisation and of "exceptional configurations" (with accidental near- zero 
modes in quenched simulations). In order to check how well some input 
operator Dq approximates chirality already, it is therefore a sensible cri- 
terion to evaluate the spectrum a{Do) and to test if it is close to this 
GW circle [99] . Let us start with the free HF in c? = 4: the spectrum of 
-Dhf in infinite volume is shown in Figure [191 cind we see that it approx- 
imates the GW circle (with p = 1) extremely well. On the other hand, 
the spectrum of the free Wilson operator covers four circles of this kind. 




(7.22) 




(7.23) 



67 



so that its real part extends up to ^ 8. (We anticipated this property 
in Section 6 when we commented on Figure [T51 ) 




O.E 1 1.E 2 2.6 3 



Re 

Figure 19: The spectrum of the truncated perfect, free HF in d = A, on 
a lattice of unit spacing and infinite size. It is very close the GW circle 
with centre and radius p = 1, which shows that it approximates chirality 
very well. 

To study this property in the presence of gauge interactions, we first 
return to the 2-fiavour Schwinger modeL We considered this model as 
described in the first part of Subsection 6.3, and we show the spectrum 
of a typical configuration at /3 = 1 / = 6 for and for Duf in Figure 
[20] [75j. In two dimensions, the spectrum of covers the vicinity of two 
circles, whereas the HF is again very close to the GW circle with p = 1. 
As an experiment, we also show in the latter case the spectrum after a 
minimal approximation to the overlap formula: we only use its first term 
in the Taylor expansion, which leads to a spectrum that can hardly be 
distinguished from the exact GW circle. As we increase the strength of 
the gauge coupling 1/a/5, the eigenvalues spread out a bit more, but 
they still follow closely the GW circle for /5 = 4 and even /3 = 2, as we 
illustrate in Figure [211 

^^In Ref. [55] we also measured the violation of the GWR directly, and we compared 
the couplings before and after the application of the overlap formula as further criteria 
for the approximate chirality of different options for Dq. 
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Figure 20: T/ie spectra of the Wilson operator (with and without a clover 
term with coefficient 1 ) and of the HF operator for a typical configuration 
in the Schwinger model at (3 = 6. The Wilson spectrum deviates strongly 
from the GW circle, whereas the HF spectrum approximates it well. In 
the latter case we also show the result if the overlap formula is approxi- 
mated by a polynomial of first order only, which is sufficient to put the 
eigenvalues quite exactly onto the GW circle fTB^ - 
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Figure 21: The spectra of the HF operator for typical configurations in 
the Schwinger model at (3 = 4 and at (3 = 2. The GW circle is still 
approximated well. We also show a polynomial correction with the Taylor 
expanded overlap formula to the first resp. second order. 

In QCD this property is more tedious to achieve. Refs. |lU7[llU2j des- 
cribe the construction of a suitable HF formulation for quenched QCD 
at P = G/g"^ = 6. It starts again from the free, truncated perfect HF and 
restores approximate criticality under gauge interaction by a link ampli- 
ficatioE0 f/j,^^ uUx,^-, u>^l. Further ingredients are a separate link 
amplification factor v for the vector term (which controls the imaginary 

^°This procedure is inspired by the method of "tadpole improvement" |108j . 
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part of the spectrum, with hardly any effect on the mass renormahsa- 
tion), and the simplest version of a "fat link"@ The latter amounts to 
the substitutions 

Ux,^l ^ (1 - + 2^ [ staples terms ] (7.24) 

for all compact link variables, where we chose a in the range 0.3. . .0.5. 
The fat link helps to pull the eigenvalues around real part 1 somewhat 
closer to the GW circle |107[ I102j . At some point, also a clover term was 
considered, but since its optimisation led to a coefficient close to we 
dropped it again. We arrived at a very satisfactory approximation for 
typical configurations at = 6, as Figure [22] shows. This plot includes 
the full spectrum on a 4"^ lattice, as well as the low eigenvalues on a 8^ 
lattice, which fill the gap near zero (this gap is generic on small lattices). 




Figure 22: The spectrum of the optimised HF operator for a typical con- 
figuration in quenched QCD at (3 = 6 on lattices of the sizes 4^ (crosses, 
full spectrum) and 8^ (squares, physical part of the spectrum). 

Later on, such a construction was also accomplished at j3 = 5.85, i.e. 
on a coarser lattice, where it is more difficult. Still this program could be 
carried out successfully |103[ 1104] . In this case, we made a compromise 
between the criteria of a minimal condition number of A' A and optimal 
locality of -DqvHF (to be discussed in this and the following Subsection). 

^^A similar parameterisation was later adapted in the approach of Refs. [lOlj . 
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An immediate consequence of the approximate chirality of Dhf is that 
the polynomial evaluation of -Dqvhf (to a fixed precision) is faster, i.e. 
the required degree is lower, as we mentioned before. In QCD we used 
Chebyshev polynomials for this purpose, which converge exponentially 
as the degree rises (see e.g. Ref. |109j )ri The required degree is then 
proportional to the square root of the condition number of the operator 
A^A. In practice one usually projects out the lowest few modes and 
treats the eigenspace spanned by them separately — this reduces the 
condition number of the remaining operator very significantly. Figure 
[23] compares these condition numbers for the Neuberger operator and 
the overlap-HF, for QCD on a 12^ lattice at /5 = 6, with A; — 1 = 1 ... 19 
modes projected out. We recognise for the overlap-HF a gain factor ^ 25, 
which we anticipated at the end of Subsection 7.3. This gain factor is 
essentially due to the reduction of the maximal A^A eigenvalue, and it 
persists practically unchanged at /5 = 5.85; details can be found in Ref. 
[lOi] (Table 1). 

7.6 Locality and rotation symmetry 

Now we move on to the point where the overlap operator is already 
evaluated to a high accuracy of at least 10~^^ (in many case, the precision 
was also set to 10~^^ or 10~^^). First we compare the level of locality 
of the overlap-HF to the standard formulation D-^. A strong gain in 
this level was first observed for the free fermion in Ref. ^99j, which was 
one of the motivations to generalise the overlap operator, as described in 
Subsection 7.3@ In Figure [2^ we show this property in c? = 2: on top 
we see that the couplings in Z^ovHF decay much faster than in D^- Since 
the -DovHF couplings follow closely a single exponential curve, this plot 
also illustrates an improved rotation symmetry. Both properties can be 
understood based on the ultralocality and the good rotation symmetry 
of -Dhf, together with relation (17.211) . (Of course, ultralocality also holds 
for but the relation (17.2 ip does not, hence this property is not 

approximately inherited in D^.) In fact, the level of locality is related 
to the minimal separation of a Dq eigenvalue from p (the centre of the 

*^The "minimax" polynomial provides a slightly better approximation with the 
same degree [109) . but the Chebyshev polynomial has the advantage that the use of 
huge coefRcients can be circumvented thanks to the Clenshaw recurrence formula. 
Regarding rational approximations, the Zolotarev polynomial is optimal in this case 
[TlOj. We add that in the Schwinger model it was possible to evaluate the overlap 
operator by diagonalising A, hence no polynomial was needed in d = 2. 

^•^Ultralocality, i.e. the limitation of the couplings to a finite range on the lattice 
(cf. footnote [T0|) . is impossible for any Ginsparg- Wilson operator in d > 1, as consid- 
erations in the free case show 
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Figure 23: The condition numbers of the overlap ingredient A^A = H^, 
where H is the Hermitian operator H = ^^A (and A = Dq — p, see eq. 



(7.20)), for the Neuberger operator and for the overlap-HF, in QCD on 
a Vl^ lattice at P = 6. k — 1 is the number of lowest modes which are 
projected out. The condition number of the remaining operator, i. e. the 
ratio Cfc = (largest eigenvalue of A^A)/{k^^ eigenvalue of AM), is about 
25 times lower for the overlap-HF 1102^ . from which we infer a gain factor 
of about 5 in the polynomial degree. The computational effort is roughly 
proportional to this degree. 

This gain factor is practically the same at (3 = 5.85 on a 12^ x 24 lattice 
JM. 



GW circle) |112j : this implies a link between the quests for locality of an 
overlap operator and approximate chirality of its kernel Dq. The optimal 
parameters have roughly the same trend for these criteria, but they are 
not identical (we mentioned in the previous Subsection that we made a 
compromise between maximal locality of -DqvHF and minimal condition 
number of A"!" A). 

In the lower plot in Figure [21] we see that the superior locality of 
-DovHF persists in the Schwinger model at /? = 6 |75j. In this case we 
measure the locality in the way suggested in Ref. [112] : we put a unit 
source rjy at one site y, and we consider all sites x separated from y by a 
taxi driver distance r = |a;^ — y^\ := \\x — y\\i. Then we identify the 
maximum of the norm HD^^yr/yll, which we denote as /(r), 

/(r) = ^:^[\\D^y{U)vy \\\x-yh=r]. (7.25) 
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The exponential decay of (/(r)) in r is a compelling criterion to demon- 
strate locality. For (3 = 6 we see in Figure [H] (below) that this decay is 
much faster for -DqvHF than for D-^. 



Neuberger: I Pu(x) I o 

1 (x) I + 

overlapSO-HF: I P|,(x)| □ 
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o Neuberger 
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Figure 24: The locality of the overlap-HF vs. the standard Neuberger 
operator in d = 2. On top we show the decay of the free couplings in 
the vector term and in the scalar term A (in the notation D^v = 
Ptilti + ^j, against the Euclidean distance \x\. We see that the overlap-HF 
couplings follow a much faster exponential decay, indicating a higher level 
of locality. Moreover, the couplings in the Neuberger operator are much 
more spread out, which reveals a better approximate rotation symmetry 
for the overlap-HF. 

Below we compare the locality in the Schwinger model at (3 = 6, measured 
according to eq. ^7.2^ (in the taxi driver metrics). We see that the 
overlap-HF is still by far more local. 



We proceed to QCD and first illustrate that we obtain again a higher 
degree of locality for the overlap-HF at /? = 6 [102J and aX (3 = 5.85 
|103i I104j . see Figure [251 On top, at /3 = 6 (which corresponds to a 
physical lattice spacing of about a ~ 0.093 fm)0 we still use the taxi 
driver metrics, but below {(3 = 5.85, corresponding to a ~ 0.123 fm) we 



■^^For the physical units in quenched QCD, we always refer to the Sommer scale 
|113j in this work. 
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switch to the Euchdean metrics. The observation that the decay for the 
overlap-HF is not only faster, but in the Euchdean metrics also smoother, 
confirms again that our overlap-HF formulation is both, more local and 
to a better approximation rotation invariant than the standard formu- 
lation. We add that the quality of rotation symmetry was also tested 
directly in the Schwinger model with the procedure shown before for the 
HP in Figure [T51 For the overlap-HF a smooth decay of the isotropic cor- 
relator C3 was found, similar to the HF, affirming an improved rotation 
symmetry [75\ . 



L= 12, B = 6 
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Figure 25: The locality of the overlap-HF vs. the standard operator in 
QCD. On top the decay is compared in the taxi driver metrics at P = 6, 
where we find a gain factor of almost 2 in the exponent of the decay 1102^ . 
The plot below refers to (3 = 5.85 in the Euclidean metrics, which also 
provides a comparison of the quality of rotation symmetry {1 04^ . 

At last we turn to strong gauge couplings, which correspond to rough 
configurations and therefore to coarse lattices. Generally, the overlap 
formula is only applicable to generate a valid lattice Dirac operator up to 
a certain coupling strength, where locality coUapsesj^ We see in Figure 



overlap-HF 



45 



For a theoretical discussion of this issue we refer to Ref. [114j . 
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f2E\ that the Neuberger operator is still local at (3 = 5.7 (corresponding 
to a ~ 0.17 fm), but at /5 = 5.6 no exponential decay can be observed 
anymore (for any parameter p). In contrast, the overlap-HF (where only 
the link amplification factor is adapted compared to the formulation at 
P = 5.85) is local in both cases, and at /3 = 5.7 the function (/(r)) still 
exhibits a remarkably fast decay. 

This superior locality is essentially due to the HF structure of the 
overlap kernel. By means of fat links alone the locality of the overlap 
operator can also be improved, but only marginally [115j (assuming the 
optimal value of p in each case) . 
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Figure 26: The locality of the overlap-HF (with a link amplification factor 
u and p = 1) vs. the Neuberger operator D-^, in QCD at strong coupling 
(in taxi driver metrics). At (3 = 5.7 Z)n (with an optimised parameter 
p = 1.8) is still local, but at (3 = 5.6 its locality — and therefore its 
validity as a lattice Dirac operator — collapses. The overlap-HF is local 
in both cases. At (3 = 5.7 its locality is still stronger than the one of 
at P = 6 and p = 1.4 (which is optimal for locality in that case 1112^ ). 
These measurements were done on a 12^ x 24 lattice, and the anisotropy 
causes the bending down at large r. 



7.7 The scaling behaviour 

Again referring to the perfect action background of the HF and to relation 
f l7.2ip . we also expect a good scaling behaviour for the overlap-HF. For 
the free overlap-HF, this is clearly confirmed by the dispersion relation, 
which we show for momenta p = {pi, 0, 0) in Figure [271 In contrast, 
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scales worse than the Wilson operator in this case@ Qualitatively 
the same behaviour is observed for massive overlap fermions [TS]. We also 
repeated the thermodynamic scaling tests described earlier (in Subsection 
4.1, before applying the overlap formula). The results in c? = 2, for 
three version of the overlap-HF, are by far improved compared to the 
standard overlap operator, see Figure [7S]. Here we incorporated 
the chemical potential also for the overlap fermions according to the 
prescription (I4.1ip : for an alternative method and first simulations, see 
Ref. [TT6]. 
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Figure 27: The dispersion relation of the free, massless overlap-HF ( scal- 
ing optimised version), compared to the continuum and to the standard 
overlap formulation Dn- These dispersions end when the argument of 
the square root becomes negative. To provide an overview we also include 
the dispersion for the kernel operators D^p and D^^. 

In the interacting case, we reconsidered the "meson" dispersion re- 
lations in the Schwinger model (cf. Subsection 6.3), this time for exact 
Ginsparg- Wilson operators. Also here we observe a scaling behaviour 
which is by far better for the overlap-HF than for the Neuberger oper- 
ator, as Figure [29] shows. Further scaling tests in the Schwinger model 
with the dynamical HF and the quenched re-weighted overlap-HF can 
be found in Ref. |117j (they were compared to the scaling with dynami- 
cal Wilson fermions, which has also been investigated in Ref. |118] ). In 
QCD, a systematic scaling test is tedious and still outstanding, but the 
toy model results summarised here raise optimism also in that respect. 

^^On the other hand, all GW fermions arc still free of 0{a) artifacts in the inter- 
acting case, in contrast to Dyj. 
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Figure 28: The thermodynamic scaling ratios pressure/(temperture) and 
pressure/(chemical potential)^ for free overlap fermions in d = 2. The 
hierarchy of the scaling behaviour is confirmed in all respects. 
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Figure 29: The mesonic dispersion relation in the Schwinger model with 
two types of overlap-HFs (open circles and squares). Both, the "pion" 
(on the left) and the "rj-particle" (on the right) display a scaling which 
is by far improved for the overlap-HFs compared to the standard overlap 
fermion (diamonds). Once more we obtain the best scaling by inserting 
the (scaling optimised) SO-HF into the overlap formula (squares) fTB^ . 

7.8 The link to domain wall fermions 

Finally we remark that the overlap fermion is equivalent to the domain 
wall fermion |119] in the limit of an infinite number of layers in a fifth 
direction0 It is then a practical issue if one tries to approximate this 
limit by a satisfactory number L^, or by a sufficiently precise polynomial 

'"'The meaning of this fifth direction is technical rather than geometrical, because 
it does not involve link variables. So that direction has no scale; what one needs is 
its thermodynamic limit oo. 
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approximation to the overlap formula. Also in the case of the domain 
wall fermions, the standard formulation inserts a Wilson kernel. Re- 
placing it by a HF kernel, as suggested in Ref. [99j, could lead to similar 
improvements as we demonstrated in this Section for the overlap fermion. 
The improved condition number of the kernel operator manifests itself 
here in a lower demand for Lg (for some required precision of chirality). 
Furthermore the gains in locality, approximate rotation symmetry and 
scaling are expected here in practically the same form as in the overlap 
case. 

8 Relating QCD Simulations to Chiral Per- 
turbation Theory 

8.1 Chiral Perturbation Theory 

We start this Section with a few general remarks on Chiral Perturbation 
Theory (xPT), an effective low energy theory, which we are going to 
relate to our QCD simulation results. 

When a continuous, global symmetry breaks spontaneously, we ob- 
tain a continuous set of degenerate vacuum states. Expanding around 
one selected vacuum, one distinguishes between excitations to higher en- 
ergy (which are identified with massive particles) and fluctuations, which 
preserve the ground state energy. The subgroups of the energy conserv- 
ing symmetry group can either transfer the selected vacuum to a different 
vacuum state, or leave it simply invariant. The number of generators re- 
lating different vacuum states corresponds — according to the Goldstone 
Theorem |120j — to the number of massless Nambu- Goldstone bosons 
(NGBs) involved. At low energy, the NGBs can be described by an ef- 
fective theory as fields in the coset space of the spontaneous symmetry 
breaking (SSB). Such effective theories still apply if we add a small ex- 
plicit symmetry breaking; we then deal with light quasi-NGBs, which 
dominate the low energy physics. The effective Lagrangian £efr con- 
tains terms of the quasi-NGB fields, which obey the original symmetry, 
as well as the (explicit) symmetry breaking terms. All these terms are 
hierarchically ordered according to some low energy counting rules for 
the momenta and the quasi-NGB masses. A simple example for such an 
expansion is outlined in Subsection 8.4.1. 

This concept is very general, but it was introduced in the framework 
of chiral symmetry breaking in QCD |121j . At zero quark masses the left- 
and right-handed spinors decouple (see eq. 03.121) ). so £qcd is invariant 
under their independent rotation. QCD is then assumed to exhibit a 
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chiral SSB of the pattern 



SU{Nf)L ® SUiNf) 



SU{Nf)L+R , 



(8.1) 



where Nf is the number of quark flavours involved. xPT is the cor- 
responding low energy effective theory [1211 \122\ 1123] . Following the 
general prescription, it deals with fields in the SSB coset space, U{x) G 



A small quark mass supplements a slight explicit symmetry breaking. 
The quasi-NGBs are then identified with the light mesons, i.e. the pions 
7T~^, 7r°, 7r~ for Nf = 2 — and for Nf = 3 also the kaons and 77-particles. 

In view of our lattice study, we have to put the system into a finite 
volume; we choose its shape as ^ = x T (T > L). We will refer to the 
formulations of xPT in two regimes, with different counting rules for the 
terms in Ccs[U]. The usual case — to be addressed in Subsection 8.2 — is 
characterised by a large volume, Lm^r ^ 1, where m^r is the pion mass, i.e. 
the lightest mass involved, which corresponds to the inverse correlation 
length. This is the p-regime, where finite size effects are suppressed, and 
one expands in the meson momenta and masses (p-expansion) |122] . 

The opposite situation, LniT, < 1, is denoted as the e-regime. In that 
setting, an expansion in the meson momenta is not straightforward, due 
to the dominant role of the zero modes. However, the functional integral 
over these modes can be performed by means of collective variables [123.J . 
There is a large gap to the higher modes, which can then be expanded 
again, along with the meson masses (e-expansion) [1231 11251 1126] . We 
will address that regime extensively in Subsection 8.4. 

In both regimes, the leading order of the effective Lagrangian (in 
Euclidean space) reads 



Throughout Section 8 we consider only the u- and (i-quark, and we as- 
sume their masses to be degenerated; we denote them by rrig. The co- 
efficients to the terms in C^s are the Low Energy Constants (LECs), 
and we recognise and S as the leading LECs. Experimentally the 
pion decay constant was measured as F^^ ~ 92.4 MeV. S is not directly 
accessible in experiments, but its value is assumed to be in the range 

^^Details about the symmetry groups involved and further aspects are reviewed 
extensively for instance in Ref. ^124j . 



SU{Nf) 




£eff[t/] = ^Tr[a^[/t9^f/]-isTr[^(f/ + f/t)] + ... , 
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(250 MeV)^ . . . (300 MeV)^. For instance, in the one flavour value 
around (270 MeV)^ was recently obtained based on a large A^^^ expansion 

m- 

The LECs are of physical importance, but they enter xPT as free pa- 
rameters. For a theoretical prediction of their values one has to return to 
the fundamental theory, which is QCD in this case. Due to the notorious 
lack of analytic tools for QCD at low energy, the evaluation of the LECs 
is a challenge for lattice simulations. 

The LECs in Nature correspond to their values in (practically) infi- 
nite volume, and the p-regime is close to this situation. However, these 
phenomenological values of the LECs can also be determined in the e- 
regime, in spite of the strong finite size effects. Actually one makes use 
exactly of the finite size effects to extract the physical LECs. Generally, 
we need a long Compton wave length for the pions, l/m^r, and in view 
of lattice simulations in the p-regime we have to use an even larger box 
length L. In this respect, it looks very attractive to work in the e-regime 
instead, where we can get away with a small volume |128j . 

However, such simulations face conceptual problems, which delayed 
their realisation until this century: first, to realise light pions the lattice 
fermion formulation should keep track of the chiral symmetry. Second, 
the e-regime has the peculiarity that the topology is important (in ac- 
cordance with the importance of the zero modes) jl29j . xPT predictions 
for expectation values often differ when restricted to distinct topological 
sectors, so it would be a drastic loss of information to sum them up. This 
distinction requires a sound definition of the topological charge on the 
lattice. In both respects, the use of Ginsparg- Wilson fermions is ideal, 
due to the specific properties explained in Subsections 7.1 and 8.3. 



8.2 Simulations in the p-regime 

Usually, for xPT in a finite spatial box L^, one expects the meson mo- 
menta p to be small, so that 

p^'^<^ A-kF^ . (8.3) 

The term 4,tiFt, takes a role analogous to Aqcd- Regarding the counting 
rules for the momenta and the pion mass, the condition L ^ 1/ allows 
for an application of the p-expansion |122j . It expands in the following 
dimensionless ratios, which are expected to be small and counted in the 
same order, 

1 p rriTr 



LFt, Aqcd Aqcd 



4) 



80 



In this Subsection we present quenched simulation results in the p-regime. 
In Ref. |1U4] we applied the overlap-HF (described in Section 7) at /3 = 
5.85 on a lattice of size 12^ x 24, which corresponds to a physical volume 
of ^ ~ (1.48 fm)^ X 2.96 fm (where we use again the Sommer scale [113] , 
cf. footnote m]). We evaluated 100 propagators for each of the bare quark 
masses 

arriq = 0.01, 0.02, 0.04, 0.06, 0.08 and 0.1 

(in physical units: 16.1 MeV . . . 161 MeV). At this point we re-introduce 
a general lattice spacing a. We will see that the smallest mass in this 
set is at the edge of the p-regime — even smaller quark masses will 
be considered in the e-regime (Subsection 8.4). Part of the observables 
presented in the current Subsection were also measured on the same 
lattice with the Neuberger operator at p = 1.6 (a preferred value at 
p = 5.85) [T30] . 

We include niq in the overlap operator fl7.12p in the usual way, 

I?ov(m,) = (l-^)Dov + m, , (8.5) 

which leaves the largest real overlap Dirac eigenvalue (2p/a) invariant. 
ruq represents the bare mass for the quark flavours u and d. 
We first evaluate the pion mass in three different ways: 

• TTT'n, pp is obtained from the decay of the pseudoscalar correlation 
function (P(x)P(O)), with P(x) = ipx'154'x- 

• ^n,AA is extracted from the decay of the axial-vector correlation 
function (A^i^x) A^i^O)) , with A4^{x) = ipxlblii^x- 

• ^n,pp-ss is obtained from the decay of the difference 

(P(x)P(O) - S{x)S{0)) , where S{x) = ^j^^^ 

is the scalar density. This subtraction is useful at small rriq, where 
configurations with zero modes ought to be strongly suppressed by 
the fermion determinant. In our quenched study, this suppression 
does not happen as it should, but the above subtraction in the ob- 
servable eliminates the zero mode contributions, which are mostly 
unphysical. 

The results in Figure [30] (on top) show that the pion masses follow 
to a good approximation the expected behaviour oc rriq. Deviations 
occur at the smallest masses, where we observe the hierarchy 

m^,pp > m^,AA > rn^,pp-ss , (8.6) 
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in agreement with Ref. p3]- This shows that the scalar density subtrac- 
tion is in fact profitable, since it suppresses the distortion of the linear 
behaviour down to the lightest pion mass in Figure [301 

m^^PP-ssiarrig = 0.01) ~ (289 ± 32) MeV . (8.7) 

That mass corresponds to a ratio L/^ ^ 2, which confirms that we are 
leaving the p-regime around this point. Based on the moderate quark 
masses in Figure [30|, we find a very small intercept in the chiral extrap- 
olation, 

m^,PP.ss{m, 0) = (-2 ± 24) MeV . (8.8) 
Due to quenching, one expects at small quark masses a logarithmic be- 
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Figure 30: On top: The pion mass evaluated from overlap-HFs in the 
p-regime in three different ways, as described in Subsection 8.2. 
Below: The pion masses fitted to the formula 1(8. 9\) . which is expected due 
to the logarithmic quenching artifacts in the absence of an additive mass 
renormalisation. 



82 



haviour of the form 

,2 



am: 



rrir, 



Ci + 6*2 InarUq + C^arUq , (Ci, C2, C3 : constants) . (8.9) 



Corresponding results are given for instance in Refs. [13H 1132]. Figure 
30] (below) shows the fits of our data to eq. (18.91) . which works best for 
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Figure 31: The axial current renormalisation constant Za evaluated from 
^PCAC according to eq. ^S.lCh) . For -DqvHF (on the left) we find Za ^ 1 
[TU^, in contrast to the result with the standard overlap operator 
U30f (on the right). The chiral extrapolations are given in eq. liS.ll]) . 

On the same lattice we also measured the p-meson mass in the p- 
regime using the overlap operators D-^ |130j and -DqvHF |104] . as well as 
the quark mass according to the PCAC relation, 



^PCAC 



E,(pt(x)P(o)) 



(8.10) 



where we used a symmetric nearest-neighbour difference for 84. It de- 
termines the axial-current renormalisation constant Za = "m-q/mpcAc- 
For the overlap-HF this constant is close to 1 [104j . see Figure [31] (on 
the left), which is favourable in view of the link to perturbation theory. 
This is in contrast to the large Za value found for the standard overlap 
operator |130j . see Figure ED (on the right). A chiral extrapolation leads 
to 



1.17(2) for DovHF , Za = 1.448(4) for 



(8.11) 



Regarding D^, consistent results were reported later in Refs. |133[ 1134] . 
and (somewhat surprisingly) at /3 = 6, p = 1.4 it even rises to Za — 1-55 
|134] . When one uses the improved Liischer-Weisz gauge action [135] , 
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the value of for is still in that range |132j . The application of fat 
links, however, helps to reduce Za [115]. 

As a further observable in the p-regime, we measured the pion decay 
constant by means of the relation 

^777 

based on P(a;)P(0), and based on P(a;)P(0) - S{x)S{^). The results for 
the operator -DqvHF |104] are given in Figure [321 In particular the value 
at aniq = 0.01 (the lightest quark mass in this plot) is significantly lower 
for the case of the scalar subtraction. Hence this subtraction pushes the 
result towards the phenomenological value. 
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Figure 32: The pion decay constant based on a matrix element evaluation 
in the p-regime — given by eq. li8.1S\) — using the overlap-HF I104f . 



However, a chiral extrapolation based on such data for F„ cannot 
be reliable. An extrapolated value of F^^ pp would come out clearly too 
large, as it is also the case for |130j . But in particular the instabihty 
of Ft, pp_ss at our lightest pion masses in the p-regime (obtained at 
arriq = 0.01 and 0.02) calls for a clarification by yet smaller quark masses. 
We did consider still much smaller values of rUg in the same volume. As 
the results for the pion masses suggest, we are thus leaving the p-regime. 
For the tiny masses aruq < 0.005 we enter in fact the e-regime, where 
observables like P^ have to be evaluated in completely different manners, 
see Subsection 8.4. 
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8.3 The distribution of topological charges 



In the next Subsection we are going to present result in the e-regime of 
QCD. In this Subsection we focus on topological charges, which play an 
essential role the e-regime [129j (as we mentioned before). 

As we already outlined when commenting on Figure 3, it is a priori 
not obvious how to introduce topological sectors on the lattice. However, 
if one deals with Ginsparg- Wilson fermions, a sound definition is given by 
adapting the Atiyah-Singer Theorem from the continuum and defining 
the topological charge of a lattice gauge configuration by the fermionic 
index v |86j, 

topological charge = i/ := — n_ , (8.13) 

where n± is the number of zero modes with positive/negative chirality. 
These numbers are unambiguously determined once a Ginsparg- Wilson 
Dirac operator is fixed (and in practice only chirality positive or chi- 
rality negative zero modes occur in one general non-trivial configura- 
tioio). However, for a given gauge configuration, the index for different 
Ginsparg- Wilson operators does not need to agree. Albeit the level of 
agreement should be high for smooth configurations, i.e. it should — and 
it doe^ — increase for rising values of /5. 

In view of the LEG evaluation in the e-regime, numerical measure- 
ments inside a specific topological sector — and a confrontation with the 
analytic predictions in this sector — are in principle sufficient. This 
requires the collection of a large number of configurations in a spe- 
cific sector. The (hyperbolic) "topology conserving gauge actions" S'^^p 
[SI [ISZl [M [M [110] are designed to facilitate this task, 

^hyP(f/^) = J Js%%e ^p(Up) < ^ (8.14) 

^ 1 -|-oo otherwise 

where Sp{Up) = S^^{Up) is the Wilson plaquette gauge action, and Up 
are the plaquette variables [28]. For e < l/[6(2 + v^)] ~ 0.049 topologi- 
cal transitions cannot occur under continuous deformations of the gauge 

^^A cancellation between zero modes of both chiralities is manifest for the free 
fermion, but in a random gauge background this coincidence has a probability of 
measure zero jl36j : it disappears under small variations in the gauge configuration 
without topological protection. 

The case of an undefined index also exists, but it has again probability of measure 
zero. In the overlap formula (|7.12p this corresponds to a zero denominator. Such 
configurations are "exceptional" in the topological sense [T7] , though not in the sense 
of accidental (near) zeros of D due to quenching, cf. Subsection 7.5. 

^°For instance, we observed at /3 = 6.15 on a IG'^ lattice that the index i' of Dn is 
very stable as p rises from 1.3 to 1.7; this changes less than 2% of the indices. 
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configuration [141] . But in practice we have to relax e to larger values 
to allow for reasonably strong fluctuations. For strong gauge couplings 
we can then arrange for a useful physical lattice spacing. Examples illus- 
trating the increased topological stability in the course of a Monte Carlo 
history are shown in Figure [33l 



|i = 6.18, 1/e = p = 0.8, 1/e= 1.25 p = 0.3, 1/e = 1 .52 




50 100 150 200 50 100 150 200 50 100 150 200 



conf. number conf. number conf. number 

Figure 33: Histories of the topological charge \ v\ on a 16^ lattice for three 
sets of parameters, which correspond approximately to the same physical 
scale, a ^ 0.08 fm /Hg/. 



Let us now return to the Wilson gauge action, which allows us to 
investigate also the statistical distribution of the topological charges. 
Again ai (3 = 5.85 on a 12^ x 24 lattice we compared the charges for 
the overlap-HF operator, and for the standard overlap operator at 
p = 1.6 [104] . As an example, the histories of about 200 indices for the 
same configurations are compared in Figure [311 Of course, these two 
types of indices are considerably correlated, but only about 40% really 
coincide. Nevertheless both follow well the expected Gaussian distribu- 




-10 ^ ^ ^ ^ 
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configuration number 

Figure 34: Index histories for -DqvHF (see Section 7) and for the standard 
overlap operator (at p = l.Q) for the same set of QCD configurations. 

tions oc exp(— const. ■ z/^), with a width ^ 3.3, see Figure [351 This width 
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Gaussian with same <v > and norm 



topological charge 




topological charge 



Figure 35: The histogram of Do^uf indices (on the left) and of indices 
(on the right), on a 12^ x 24 lattice at (3 = 5.85. In both cases, 1013 
configurations are included \TU^. 



fixes the topological susceptibility (cf. eq. (12.151) ) 

Xtop = :^(^') , (8.15) 

which is of importance to explain the heavy mass of the ?7'-meson |142j . 

In Figure [36] we present our results with Z^ovHF and Dn on the lattice 
referred to so far, plus a result for Dn at = 6 in the same physical 
volume (lattice size 16^ x 32). We also mark the continuum extrapolation 
according to Ref. [143] , which is fully consistent with our results. That 
measurement of Xtop was based on indices on lattices]^ The 
resulting value for xtop is compatible with the Witten-Veneziano scenario 
that much of the rj' mass is generated by a U{1) anomaly jl42j . 



8.4 The e-regime 

We first repeat that the e-regime of QCD is characterised by a relatively 
small volume, i.e. the correlation length ^ exceeds the linear box size L. 
Together with the requirement (18.31) this amounts to the conditior0 

— >L>-i-. (8.16) 

In such a box, the p-expansion of xPT fails, in particular because of 
the importance of the zero modes. However, the latter can be treated 
separately by means of collective variables, and the higher modes — along 



A compilation of earlier lattice results for xtop (with various methods) is given in 
Ref. [14] . 

^^In such small volumes the notion of ttItt is problematic, but this inequality can be 
interpreted by referring to the would-be pion mass in an extended volume. 
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Figure 36: The topological susceptibility measured by indices of Do^uf 
and of in o volume V = (1.48 fm)^ x 2.96 fm, with two different 
lattice spacings a. Our data — given in detail in Ref. [104^ — agree well 
with the continuum extrapolation in Ref. 



with the pion mass — are then captured by the e-expansion |123] . One 
now counts the ratios 



p" 



1 



A 



QCD 



^^QCD 



(8.17) 



as small quantities in the same order. 

This setting, where pions are squeezed into a tiny box, cannot be 
considered physical. Nevertheless there is a strong motivation for its 
numerical study: the point is that the finite size effects are parametrised 
by the LECs of the effective chiral Lagrangian as they occur in infinite 
volume. Hence the physical values of the LECs can (in principle) be 
evaluated even in such an unphysically small box, as we mentioned in 
the introduction to Section 8. 

Unfortunately, most results for the LECs in the e-regime are obtained 
in the quenched approximation so far, hence they are affected by (mostly 
logarithmic) finite size effects |144] . So the final results by this method 
still have to wait for the feasibility of QCD simulations with dynamical, 
chiral quarks. These prospects will be commented on in Subsection 9.2. 



8.4.1 A 3-loop calculation in the e-expansion 

The chiral symmetry breaking for two fiavours is locally isomorphic to 
orthogonal groups. 



SU{2)®SU{2) SU{2) 



0(4) ^ 0(3) . 



(8.18) 
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Generally, the spontaneous symmetry breaking 0{N) 0{N — 1) gener- 
ates N—1 NGBslffl They can be described by the non-linear cr-model and 
studied with the formalism of xPT. In d = 3 and d = 4 the corresponding 
Lagrangian involves an infinite string of terms, which can be ordered by 
the power of the momenta involved. In a volume V = L'^ (with periodic 
boundary conditions) we count the momenta as (9^ = 0(L"^), and we 
obtain the leading terms 

+ ^gf\d,SdJf + 0{L-') , S\x)^l, (8.19) 

where we attach an independent LEG to each term. We may also add 
terms that break the 0{N) symmetry explicitly through a small, constant 
magnetic field H, which adopts the role of the light quark masses, 

+ h^^l{HS){d^Sd^S) + hfl{Hd^S) + ... . (8.20) 

Then we assemble the total Lagrangian 

We count H := \H\ = 0{V-^), so that the NGBs with = EH/F^ 
feel the finite size strongly, rriL <^ 1. Then the partition function can be 
e-expanded in the dimensionless ratio e = L^^'^/F^ in = 3 and 4. 

To compute the free energy = — In Z one first has to find a way to 
handle the zero modes. We mentioned before that this can be achieved 
by using collective variable (123] . In the present case, one considers the 
magnetisation 

M := J d'^xS{x) = \M\ne , (8.21) 

where e = H/\H\ is a fixed unit vector and Q G 0{N) is integrated 
over in the functional integral. (The corresponding unitary integration 
for quark flavours is discussed in Refs. [145] .) Here the vector field 
S can be decomposed into one dominant component plus fluctuations 
7r(x) = (7ri(x), . . . 7rAr_.i(a;)), TTi{x) = 0{L^~'^^'^), which are treated per- 
turbatively. In this way, the partition function can be evaluated order 
by order. 

This was first carried out to 2 loops in Ref. jl26j , where the func- 
tional measure was treated with the Faddeev-Popov procedure (the con- 
straint (5(5"^ — 1) is implemented in an exponential form). A different 

^•^There is, however, no other symmetry breaking 0{N) 0{n), N > n, which is 
locally isomorphic to SU{Nf) (g) SU{Nf) SU{Nf) for any Nf. 
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method was applied in Ref . |146] , which was based on the Polyakov mea- 
sure |147] . That measure was originally introduced in the framework of 
string theory and it captures the above constraint by integrating T>n over 
the e-expanded elements of a metrics in flavour space. In this way we 
computed the partition function to 3 loops. In both cases dimensional 
regularisation [31] was used. It turned out that the free energy is in 
fact perturbatively renormalisable order by order (although the number 
of required counter terms increases rapidly in each order). This property 
is highly non-trivial due to the requirement that the renormalised LECs 
must not pick up any volume dependence. This is realised both, in d = 3 
and in d = 4, for the large number of terms occurring to the 3-loop or- 
der, thanks to numerous cancellations between forbidden contributions. 
These calculations represent therefore a sensitive test for the validity of 
the e-expansion scheme, as well as the methods used for the treatment 
of the functional measure. 

Before chiral lattice fermions became available, the program of evalu- 
ating LECs through simulations in the e-regime was tested in the frame- 
work of this spin model [148] . 

8.4.2 The chiral condensate 

Chiral Random Matrix Theory (RMT) conjectures predictions for the 
low lying eigenvalues, ordered as A^, n = 1, 2, 3 . . . (excluding possible 
zero eigenvalues) of the Dirac operator in the e-regime; for a review, see 
Ref. |149j . More precisely, it conjectures densities of the dimensionless 
variables HVXn, where S is the chiral condensate in the effective La- 
grangian (18.21) . Here we focus on the variable z := TiVXi^p, where Ai,p 
emerges from the leading non-zero eigenvalue Ai if the spectral circle 
of the overlap operator is mapped stereographically onto the imaginary 
axis, Ai,p = |Ai/(l - Ai/2p)|. 

These RMT predictions depend on |z/|, the absolute value of the topo- 
logical charge. Here we make use of the explicit formulae [150] for the 
density of the first non-zero (re-scaled) eigenvalues z in the sectors 
which we denote by p^}'^^\z). For the lowest eigenvalues, the particular 
density p^^^ was first confirmed by staggered fermion simulations (results 
are summarised in Ref. |149j ). But in those studies the charged sector^ 
yielded the very same density, in contradiction to RMT. 

The distinction between the topological sectors was first observed to 
hold for Djss, following the RMT predictions to a good precision [152] . 
if the linear box size exceeds a lower limit of about L>1.1 fm (the 

^ In that topological charge has to be introduced by some traditional method 

like cooling [151] , 
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exact limit depends on the criterion, of course) @ The predictions for 
the densities of the leading non-zero eigenvalues in the sectors |z/| = 0, 1 
and 2 are shown in Figure [371 on the left; we see that zero modes repel 
the finite eigenvalues. On the right-hand-side of Figure [37] we present 
results for the corresponding cumulative density with on a 10'' lattice 
at (3 = 5.85 (box length ^ 1.23 fm). Once the predicted densities Pi'^'^ 




Figure 37: On the left: RMT predictions for the leading non-zero Dirac 
eigenvalue in the topological sectors with charge |z/| = 0, 1 and 2. 
On the right: RMT predictions (lines) and simulation results for the cor- 
responding cumulative densities. The numerical data were obtained with 
Dt<s on a 10"^ lattice at [3 = 5.85 (see first work in Ref. flS^ ). and they 
do essentially follow the RMT predictions. The confidence level for the 
agreement of cumulative densities has been verified with the Kolmogorov- 
Smirnov test 1109^ . 

(marked by lines) are well reproduced, we can read off the value of S, 
which is the only free fitting parameter for all topological sectors. We 
proceed to larger lattices and show our result for Z^ovHF on a 12'^ x 24 
lattice at /5 = 5.85 in Figure [38] |104] . The optimal fit shown in this plot, 
and its counterpart for the Neuberger operator yield 

S^/^ = 298(4) MeV (from DovHf) , ^^^^ = 301(4) MeV (from D^) . 

(8.22) 

These fits focus on the lowest eigenvalues resp. energies, where chiral 
RMT is most reliable. Clearly, in this range the neutral sector (z/ = 0) 
dominates. In the case of -Dn , the charged sectors |z/| = 1 and 2 alone 
would favour a different S value |104j . Such ambiguities are quite strong 
in the results with smeared staggered fermions [153] . In the case of -DqvHF, 
however, a unique S works well for all the three sectors |z/| =0, 1, 2, up 
to about Z ~ 3, clS Figure [38] shows. This range extends well beyond 

^^Meanwhile a topological splitting has also been observed to set in for staggered 
fermions if the link variables are strongly smeared |153j . 
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the Thouless value ^Thouiess = F^W < 1, which is often understood as a 
threshold for the RMT applicability. 



12^x24, p=5.85 




z 



Figure 38: The cumulative density of the (Mohius projected) lowest Dirac 
eigenvalue Xi^p of the overlap-HF operator, in the topological sectors |z^| = 
0, 1 and 2. We compare the chiral RMT predictions to our data for 
z = "EVXi^p with T,^^^ = 298 MeV — the optimal value in the neutral 
sector (u = 0). This value works well up to z ^ 3 for all topological 
sectors, i.e. well beyond the Thouless value -Zxhouiess ^ 1; which is often 
considered a theoretical bound for the applicability of these predictions. 

As an alternative approach to test the agreement of our data with 
the chiral RMT, and to extract a value for E, we now consider the mean 
values of the leading non-zero Dirac eigenvalues Ai in all the charge sec- 
tors up to |z/| = 5. In physical units, the results (Xi^p) agree very well 
for different overlap operators and lattice spacings — see Figure [39] — 
although this consideration extends beyond very low energy. Each single 
result for (Ai^p)|i,| can then be matched to the RMT value for a specific 
choice of S. It is very remarkable that all these 18 results are compatible 
with RMT if we choose 

S = (290(6) MeV)=^ , (8.23) 

as Figure 15^ also shows. 

A renormalisation procedure for S obtained in this way is discussed 
in Ref. |154j . However, we will only use this quenched lattice result as a 
fitting input in Subsection 8.4.3, so we stay with the bare condensate S 
for our fixed lattice parameters. 

We add that there are also first applications of this technique with dy- 
namical overlap fermions [155\ . 
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Figure 39: The mean values of the first non-zero Dirac eigenvalue (in 
physical units) in the charge sectors |z/| = 0...5. All these numerical 
results 11U4]I agree with chiral RMT if we choose T}'^ = 290(6) MeV. 



8.4.3 The pion decay constant determined from the axial cur- 
rent correlator 

In order to relate our quenched simulation results to the effective low en- 
ergy theory, we now refer to quenched xPT. In that framework, mesonic 
correlation functions were calculated to the first order in Refs. |156[ 1157] . 
The vector current correlation function vanishes, while the scalar and 
pseudoscalar correlators involve already in the first order additional, 
quenching specific LECs, which obstruct the access to the physical LECs 
in the Lagrangian (18.21) . Therefore we first focus on the axial- vector cor- 
relator, which only depends on E and in the first order. In particular 
we are going to compare our data to the quenched xPT prediction in 
a (periodic) volume V = x T |157] . for the topological sectors with 
charge ±z/. 



Zl-{A,{t) A4(0))h 
hi{r) 



2( ^ + 2m,S|,|(^,)T/ii(r; 



r — r 



T 



Iu{Zq)Ky{Zq) + Iy+i{Zq)K^^i{Zq 



(8.24) 



V 

+ — 



where 



Ai{t) = ^ ^(t, x)7574^/'(t, x) 



(t>0) 



(8.25) 



93 



is the bare axial-vector current at 3-momentum p = 0. I^, and are 
modified Bessel functions, and Zg := HVrUg (in analogy to the dimen- 
sionless variable z in Subsection 8.4.2). 

It is remarkable that this prediction in the e-regime has the shape of 
a parabola with a minimum at t = T/2. This is in qualitative contrast to 
the cosh behaviour, which is standard in large volumes. S affects both, 
the curvature and the minimum of this parabola, whereas only appears 
in the additive constant — that feature is helpful for its evaluation. 

A first comparison of this curve to lattice data was presented in Ref. 
[158] ■ using at /3 = 6, p = 1.4, arrig = 0.01 on lattice sizes 10^ x 24 and 
12^. For the anisotropic volume the linear size of L ~ 0.93 fm turned out 
to be too small: the data for (A^lt) A4^{0)) 1^2 were incompatible with the 
parabola of eq. (18.241) for any positive S. This observation is consistent 
with the lower bound for L that we also found for the agreement of 
the microscopic spectrum with chiral RMT (L>1.1 fm, see Subsection 
8.4.3). 

Another observation in that study was that the corresponding Monte 
Carlo history in z/ = is plagued by strong spikes, giving rise to large 
statistical errors. A huge statistics of 0(10"^) topologically neutral config- 
urations would be required for conclusive results. These spikes occur for 
the configurations with a tiny (non-zero) Dirac eigenvalue Xi^p. It agrees 
again with chiral RMT that such configurations are most frequent in the 
topologically neutral sector, see Figure ETl (A method called "low mode 
averaging" was designed and applied to alleviate this problem |159j .) 

However, without applying that method we obtained a decent agree- 
ment with the prediction (18.241) in the 12^ lattice mentioned above {V ~ 
(1.12 fm)^) in the sector |z/| = 1 [158J. In view of the leading LECs, it 
seems unfortunately impossible to extract a value of E from such data 
(although it is encoded in Zg), since the theoretical curvature depends 
on it only in an extremely weak way@ On the other hand, can be 
extracted quite well from the vicinity of the minimum at t = T/2, but 
the value found in Ref. [158] was too large. 

Next a study of that kind appeared in Ref. |133j . which also used Dn, 
at /5 = 5.85 (and p = 1.6), now on a 10^ x 20 lattice. This work analysed 
the sectors =0 and 1. As a reason for this limitation the authors 
referred to the condition |z/| ^ (^^)- As we mentioned in Subsection 
8.4.2, one expects (z/^) oc V (up to artifacts), see eq. (18.151) . hence this 
limitation was imposed by the volume. 

In Ref. |104] we presented again results at /3 = 5.85 on a 12^ x 24 
lattice, where the volume admits =2. We measured for both, -DqvHF 

^®Only in the sector v — the sensitivity to E is significant, but there we run into 
the statistical problem mentioned before. 
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arrig 


-DovHF 




0.001 


= (110 ± 8) MeV 


F^ = (109 ± 11) MeV 


0.003 


= (113 ±7) MeV 


F^ = (110 ± 11) MeV 


0.005 


F^ = (115 ± 6) MeV 


F^ = (111 ±4) MeV 



Table 2: Our results for the pion decay constant F^^, evaluated in the 
e-regime based on the axial-current correlation function \8.2Ji\j . These 
results are obtained at f5 = 5.85 on a 12^ x 24 lattice. The statistics 
involved 100 propagators for at arrig = 0.005, and 50 propagators in 
all other cases. F^^ was determined from fits to the quenched xPT formula 
1^8.24^ in the range t/a& [11,13]. 



and Dn, the axial-vector correlators at the masses aruq = 0.001, 0.003 
and 0.005, which are safely in the e-regime. We fitted the data to eq. 
f l8.24p by inserting the chirally extrapolated factors Za (1.17 for Z^ovHF 
[104J and 1.45 for [130] . see Subsection 8.2), along with the S values 
in eq. (18.221) . For each of the overlap operators we performed — at each 
of the quark masses — a global fit over the topological sectors that we 
considered. The result for -DqvHF is shown in Figure HOI It revealed for 
the first time a convincing distinction between the sectors |z/| = 1 and 
|z/| = 2 — this predicted topological splitting could not be observed for 
Z^N up to now. For at anig = 0.005 we also include the neutral sector; 
as expected it has clearly larger errors than the charged sectors, but it 
is helpful nevertheless to reduce the error on Fj^ in the global fit. The 
values for F^, obtained in this way for Z^ovHF and for are in accurate 
agreement, as Table [2] shows. 



8.4.4 The pion decay constant based on the zero modes 

At last we review our results based on an alternative method to evaluate 
in the e-regime |104] . This method was introduced in Ref. [160] . 
and it involves solely the zero mode contributions to the pseudoscalar 
correlation function. Here one works directly in the chiral limit. Let us 
briefly summarise the main idea of this approach. 

Ref. [160j computed the chiral Lagrangian to the next-to-next-to- 
leading order in quenched xPT, -^^q^^pT- It can be written in a form 
that involves an auxiliary scalar field $0; which is coupled to the quasi 
Nambu-Goldstone field U by a new LEG denoted as K. The auxiliary 
field also contributes 

^^'^[*o] = T^W,$o + #$^ (8.26) 
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Figure 40: Lattice data obtained with -DqvHF for the axial-current corre- 
lation functions, measured separately in the topological sectors |z/| = 1 
and 2. The curves are global fits (over both sectors) to the quenched xPT 



formula ( 8.24 h for each of our three masses in the e-regime. They single 
out the values of Ft, given in Table\^ 



to -^^q^PT) which brings in ao and mo as another two quenching specific 
LECs, in addition to K. The field $o supplements the quenching effects; 
in the dynamical case it decouples form the field U . 

It is ambiguous how to count these additional terms in the quenched 
e-expansion. Ref. [160] assumes the action terms with the coefficients 
tto and K\/Wc to be of 0(1), whereas the one with mo is in 0(e). The 
last assumption is somewhat unusual (for instance, it differs from the 
framework of Subsection 8.5.1). Nevertheless it is an acceptable possi- 
bility, which simplifies this approach since it removes the auxiliary mass 
term from the dominant order. If one further defines the dimensionless 
parameter 

AN'^KF 

(8.27) 



a 



then only the LECs and a occur in this order. 

For Nf valence quark fiavours, this approach considers the correlation 
function of the pseudoscalar density P{x) (defined in Subsection 8.2), 
which is decomposed into a connected plus a disconnected part. In a 
spectral decomposition of the propagators one obtains the residuum in 
terms of the zero modes, 



lim 



{m,Vf{P{x)Pm. 



connected : C'^l^^{x) 
disconnected : 0|^| (x) 



NfC[^^{x) + N}c[^^{x) 
-{v]{x)vu{x)-vl{Q)v,m\,\ 
{v]{x)v,{x)-vl{Q)vum\.\ ■ (8.2^ 
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The vectors Vj denote the (exact) zero modes of the Ginsparg- Wilson 
operator at rriq = 0. In the terms for C^^^ these zero modes are summed 
over. 

Next we consider the spatial integral J d^x P{x)P{0). Now the above 
procedure for the correlation function leads to functions C|^*j(t), i = 1, 2, 
which are given explicitly in Ref. |16U] . In principle, these functions could 
be measured and fitted to the predictions in order to determine and 
a. In practice, however, it is much more promising to consider just the 
leading term in the expansion at t = T/2, 

^,jClp)\,=r/2 = Di;^s + 0{s') , s = t-^, .= 1,2. (8.29) 

The explicit slope functions D^^J^ in a volume V = x T are given in 
Refs. |160l 11031 1104j . (They also involve a shape coefficient, which we 
computed for our anisotropic volume according to the prescription in Ref. 

m-) 

We evaluated the LECs and a from fits to the linear term in 
eq. fl8.29p . For each of our lattice sizes and for each type of overlap 
operator we performed a global fit over the topological sectors |z/| = 1 
and 2, in a fitting range Smax • The slopes tend to be stable over a variety 
of fitting ranges s G [— Smax? Smax], Smax = CL, 2a, 3a.... The deduced 
optimal values for are shown in Figures HI], and the values for F-n- 
and a at s^g^^/a = 1 are given in Table [31 We see that the results 
for different lattice spacings and overlap Dirac operators are in good 
agreement. Considering also a(sinax), we found the most stable plateau 

for DovHF HDl. 



Dirac operator 


DovUF 








5.85 


5.85 


6 


lattice size 


12^ X 24 


12^ X 24 


16^ X 32 


F. 


(80 ± 14) MeV 


(74± 11) MeV 


(75 ± 24) MeV 


a 


-17± 10 


-19 ±8 


-21 ± 15 



Table 3: Our results in the e-regime for the pion decay constant F^^ — 
along with the quenching specific LEC a given in eq. 1^8.27^ — based on 
the zero mode contributions to the pseudoscalar correlation function. We 
give results for the fitting range Smax = ci, which is most adequate in the 
light of eq. ^KM) . 

The values that we now obtain for F^ is below those of Section 8.4.3, 
which used a different observable and a different e-counting rule for the 



97 



200 p 

180- 

160- 

„ 140- 
> 

" 120- 
100- 
80- 
60- 

4oL 



12''x 24 p=5.85 overlap-HF 
12^x 24 p=5.85 standard overlap - 
1 6^x 32 p=6 standard overlap 



2 

s / a 



Figure 41: The results for F^^ based on a global fit of our data to a 
quenched xPT prediction for the zero mode contributions to the pseu- 
doscalar correlations function. We show the results of a two parameter 
fit (for Ft, and a) over the ranges s G [T/2 — Sma^,T/2 + Smax]- 



quenched terms. In fact, the results in Tableware close to the phe- 
nomenological value (the latter is shifted down to ~ 86 MeV if one ex- 
trapolates to the chiral limit |161] ). Our values for and a obtained 
from the zero modes are somewhat below the values reported in Ref. |160] 
based on the same method. We suspect that the anisotropic shape of our 
volumes, T = 2L, could be the main source of this deviation |162j . 



9 Epilogue 

The central concept in this work are block variable RGTs (renormalisa- 
tion group transformations) applied to lattice regularised quantum field 
theories. We described this method for various types of field theoretic 
models, and focused in particular on the limit obtained under iterated 
RGTs. This leads to lattice formulations with fascinating properties: in 
particular the symmetries and the scaling quantities of the continuum 
theories can — in principle — be reproduced exactly on the lattice. This 
amazing feature also includes exact supersymmetry and precise topo- 
logical sectors on the lattice. The corresponding perfect lattice actions 
were constructed explicitly for the case of free particle. For interacting 
fields this can in general only be achieved approximately (as an excep- 
tion we considered the Gross- Neveu model in the large N limit). We 
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discussed perturbatively perfect actions in various models, as well as is- 
sues of truncation, parameterisation and gauging. We also demonstrated 
that different kinds of anomalies are represented correctly in this way. 

We summarised a variety of simulation results based on such ap- 
proximately perfect actions. They reveal a number of properties which 
are superior to the standard lattice formulations, in particular an im- 
proved symmetry and scaling behaviour. The chiral symmetry of mass- 
less fermions can be rendered exact by means of the overlap formula. This 
provides a lattice formulation, which performs very well in toy model sim- 
ulations, and it is currently being applied to QCD. 

9.1 The status of overlap-HF applications in QCD 
with hght quarks 

In Section 7 we reviewed our construction of overlap hypercube Dirac 
operators -DqvHF, which are especially suitable at lattice spacings of a ~ 
0.093 fm and a ~ 0.123 fm. In both cases, they display a strongly 
improved locality compared to the standard overlap operator D^. Hence 
-DqvHF defines chiral fermions on coarser lattices. 

Section 8 summarised quenched simulations with -DqvHF and with 
in a volume V ~ (1.48 fm)^ x (2.96 fm) at /3 = 5.85 and at /3 = 6. 

Subsection 8.2 dealt with the p-regime, where we measured the meson 
masses and mp (though only the former was presented here, for ijip 
we refer to Ref. [104] ). the quark mass mpcAc (based on the axial Ward 
identity) and the pion decay constant F^^ at bare quark masses ranging 
from 16.1 MeV to 161 MeV. The results for m^r and nip are similar for 
-DqvHF and D^. On the other hand, for -Dqvhf the quark mass mpcAc is 
much closer to iriq than in the standard overlap formulation. This implies 
Za ~ 1, which is favourable for the connection to perturbation theory. 
Regarding F^, it turned out that the data obtained in the p-regime can 
hardly be extrapolated to the chiral limit. 

In Subsection 8.3 we discussed a way to stabilise the topological sec- 
tor in a Monte Carlo history. We also considered the distribution of 
a large number of topological charges defined by the fermion indices of 
-DqvHF or of Dn- We found histograms which approximate well a Gaus- 
sian distribution, consistent with the conservation of parity invariance. 
The resulting topological susceptibility is in good agreement with the 
literature, and with the Witten-Veneziano scenario. 

In Subsections 8.4 we proceeded to the e-regime, where we first sum- 
marised a 3-loop calculation in the framework of the e-expansion. Nu- 
merically we determined a value for the chiral condensate from the distri- 
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bution of the lowest eigenvalues. For both, Z^ovHF and we identified 
S close to (300 MeV)^. 

We evaluated F^^ in the e-regime in two ways, from the axial-current 
correlation and from the zero mode contributions to the correlation of 
the pseudoscalar density. These two methods handle the e-counting of 
the quenched terms differently, and they yield different values for F^^. 
The axial-current method leads to F^^ ^ 110 MeV, which is consis- 
tent with part of the quenched results in the literature. The zero mode 
method leads to a lower Fj^, in the vicinity of the phenomenological value. 
The result of Ref. |159j — empolying yet another method, based on the 
A J =1/2 rule, still in the e-regime — is in between. 

From the current results we conclude that the methods applied here 
do have the potential to evaluate at least the leading LECs from lattice 
simulations in the e-regime. The quenched data match the analytical 
predictions qualitatively (if the volume is not too small) and — in the 
setting we considered — they lead to results in the magnitude of the 
LECs in Nature. However, the quenched results are volume dependent 
and in addition ambiguous: different methods yield different values. 

For values that can be confronted with phenomenology in detail, 
simulations with dynamical quarks will be needed. In particular the e- 
regime requires then dynamical Ginsparg- Wilson fermions. This regime 
is promising in view of the lattice size. Also the option to extract physical 
information from single topological sectors is attractive, since it is very 
difficult to change the sector frequently in the course of Hybrid Monte 
Carlo histories. The question is if one is able to handle sufficiently small 
quark masses in dynamical simulations. 

9.2 Prospects for dynamical simulations with chiral 
fermions 

We just pointed out that the results for light quarks can be linked suc- 
cessfully to xPT (Chiral Perturbation Theory), but to some extent they 
are obstructed up to now by the quenched approximation. Quenching 
has been necessary so far in QCD simulations with chiral quarks due to 
limitations in the computational resources. 

We hope to overcome this limitation — and the systematic errors 
that it causes — in the foreseeable future, and to be able to proceed 
to simulations with dynamical chiral quarks. However, in addition to 
powerful machines this step also requires new algorithmic tools, which 
are currently under consideration. Also in this respect overlap-HFs open 
up new perspectives. 
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Due to the similarity with the hypercubic kernel, a low polynomial 
of this kernel can be used as a numerically cheap way to evaluate the 
fermionic force in the Hybrid Monte Carlo (HMC, [78j) algorithm. This 
approach is conceptually guided by the algorithm that we used for the 
simulation of quasi-perfect staggered fermions (see Ref. [77] and Subsec- 
tion 6.3.2). Also for the overlap-HF formulations of Ref. [7S] it is under 



investigation in the two-flavour Schwinger model. First the extreme case 
of using directly D^f in the force term was studied in Ref. [163] . which 
reported a decreasing acceptance rate for an increasing volume. However, 
it turned out that the acceptance rate can be improved by an order of 
magnitude by correcting the fermionic force term at least to a low accu- 
racy like 0.005, which allows for efficient HMC simulations with very light 
(degenerate) overlap-HFs [164j . Of course, the overlap operator has be to 
very precise in the Metropohs accept/reject step (we set it to 10~^^). We 
verified the acceptance rate as well as algorithmic requirements, namely 
area conservation and reversibility. 

In our simulations at /? = 5 on a 16 x 16 lattice [164] (with the 
plaquette gauge action) we confirmed the similarity between the kernel 
and the overlap operator by considering the spectra, and we observed 
a high level of locality for this formulation. The characteristic decay 
hardly changes in the fermion mass range 0.03 . . . 0.24, and the locality 
is by far superior to the standard overlap operator. As an observable 
we evaluate the chiral condensate at light fermion masses, based on the 
ratio between low lying Dirac eigenvalues in different topological sectors 
(the relevant formulae were taken from Ref. pTOS^J). This represents one 
of the first measurements with dynamical overlap fermions. Our results 
for the chiral condensate at various masses are in excellent agreement 
with analytic predictions, which were obtained with bosonisation and 
low energy approximations 



In QCD, simulations with dynamical overlap fermions are still in an 
early stage; for recent status reports we refer to Refs. [167] ^ Regarding 
the approximately chiral fermion formulations explained in this work, 
dynamical HFs are currently investigated, and interesting results for their 
phase diagram are available already [169] . It provides access to a mass 
ratio m^/mp < 0.8 at the thermal crossover, which is not the case for 
Wilson fermions. In this context we add that the truncated classically 
perfect action of Ref. [13], and the "chirally improved" formulation of 
Refs. |101j . are currently applied in dynamical simulations as well |170] . 



^^After the first version of this work was written, a remarkable study appeared 
in Ref. [168] . It evahiates the chiral condensate from simulations with dynamical 
overlap fermions in the topologically trivial sector, applying the technique described 
in Subsection 8.4.2. 
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Also the topology conserving gauge actions — that we discussed in 
Subsection 8.3 — appears to be helpful in simulations with dynamical 
quarks due to the suppression of small plaquette value, as suggested in 
Ref. p8] and tested in Ref. [TiO] . 

The quenched studies reported here show that we have methods at 
hand, which are applicable for instance for the evaluation of Low Energy 
Constants in the chiral Lagrangian from first principles. These constants 
play a prominent role in QCD at low energy, hence their determination 
is a major challenge and a sensitive test for QCD, to be addressed with 
dynamical Ginsparg- Wilson fermions. 

10 Further Fields of Research 

In my diploma thesis I investigated mathematical and mechanical models 
for certain brain activities and the resulting neural signals. More pre- 
cisely, I analysed numerically and analytically the time evolution of jerky 
eye motions and the required coordination of nerve signals, particularly 
in view of the non-Abelian dynamics. 

For short periods I was involved in two experiments at CERN that 
observed the electroweak mixing angle (CHARM II) and neutral kaon 
decays (CP-LEAR). 

In a work with J.J. Giambiagi we discussed solutions of the spherically 
symmetric wave equation and Klein-Gordon equation in an arbitrary real 
number of spatial and temporal dimensions [171] . Starting from a given 
solution, we presented various procedures to generate further solutions 
in the same or in different dimensions. The transition from odd to even 
or non-integer dimensions could be performed by fractional derivation or 
integration. We also discussed the analytic continuation to arbitrary real 
powers of the D'Alembert operator. Finally we worked out operators 
which transform a time into a space coordinate and vice versa, and we 
commented on a possible apphcation to black holes. 

In a completely different project we addressed statistical mechanics 
and investigated the phase diagram of the 2d 0(3) lattice model with a 9- 
term [172j . The simulation made use of the Wolff cluster algorithm |173] 
and required the construction of an improved estimator, which could only 
be achieved by a new method. The clusters carry half integer topological 
charge and most of the charged clusters can be identified with merons. 
We confirmed to a high accuracy the prediction by I. Affleck et al. [174] , 
who expected a second order phase transition at 6* = vr due to analogies to 
Haldane's conjecture about quantum spin chains [175] . The traditional 
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picture explains this phase transition by the neutrahsation of vortex- 
antivortex pairs. In the framework of our "meron cluster algorithm", 
this picture could be formulated precisely, without breaking the 0(3) 
symmetry. 

Further algorithmic elements were added in Ref. [176 j. Recently an- 
other meron cluster simulation was performed successfully for the XY 
spin chain [TB] . 

For a few years, I have been working in a collaboration which studied 
large N reduced matrix models as candidates for a non-perturbative for- 
mulation of string theory. In particular, we have simulated variants of the 
so-called IKKT model (or IIB matrix model) |177j . which corresponds to 
a dimensional reduction of super Yang-Mills theory to zero dimensions. 
This model is formulated in 10 (target) dimensions and it has a complex 
action, which can hardly be simulated. We studied its 4d counterpart 
[178] . which has a real positive action, as well as the lOd version with 
a quenched phase [179] (a summary is given in Ref. |180] ). The space- 
time coordinates arise dynamically, and we found numerically that their 
distribution does not support the scenario of spontaneous Lorentz sym- 
metry breaking, which could single out a subset of extended dimensions. 
(This question was also addressed in Refs. |181j .) We did see, however, a 
well-defined large N scaling, and we observed the area law for the Wilson 
loop to hold in some regime. 

In the context of large reduced matrix models, we are now study- 
ing field theory on spaces with a non- commutative (NC) geometry. This 
was the subject of the Ph.D. thesis of my collaborators A. Bigarini, F. 
Hofheinz, Y. Susaki and J. Volkholz. A (fuzzy) lattice version of 2d U{1) 
gauge theory on a NC plane can be mapped onto a twisted Eguchi-Kawai 
model [18211183] . This mathematically rigorous map enables Monte Carlo 
simulations. We observed numerically the large A^ "double scaling" of 1- 
and 2-point functions of the Wilson loops. This corresponds to a simul- 
taneous continuum and infinite volume limit at fixed non-commutativity. 
Thus we obtained the first fully non-perturbative evidence in favour of 
the renormalisability of a NC field theory [18411185] . The area law for the 
Wilson loop holds at small area as in the commutative U{N —>■ oo) model 
[186] , but at large area we observed instead a linearly rising complex 
phase. This is consistent with an Aharonov-Bohm type picture known 
from the Seiberg-Witten map to string theory, and from 2d models for 
the quantum Hall effect. The difference is, however, that we discovered 
this feature dynamical effect. 

For pure U{N) gauge theories in a commutative plane, the expec- 
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tation values of Wilson loops only depend on their (oriented) area, i.e. 
they are invariant under area preserving diffeomorphisms. This prop- 
erty enables analytic solutions. Hence the question arises if the same 
holds for NC 2d gauge theories. First perturbative calculations seemed 
to confirm this |187] . but proceeding to higher accuracy a reduction of 
the symmetry group to SL{2,R) was observed [188] . However, detailed 
numerical measurements show that non-perturbatively this symmetry is 
broken completely, including the SL{2,R) subgroup |189j . Hence ana- 
lytic solutions are unlikely in that case, but the theories may display a 
rich structure. 

Next our numerical pilot studies of NC field theory addressed the Xcj)"^ 
model in two and three dimensions |19m I185[ 1191] . The results reveal 
an interesting phase diagram, and in particular the occurrence of a new 
"striped phase", as it was first conjectured in by Gubser and Sondhi 
|192] . Its existence shows that the famous UV/IR mixing of divergences 
in NC field theory is not an artifact of perturbation theory, but it belongs 
to the very nature of the system. In three dimensions (with a commuta- 
tive Euclidean time) we could again safely extrapolate our results to the 
double scaling limit. In particular the dispersion relation stabilises, and 
it reveals an IR divergence in part of the phase diagram. Although this 
striped phase involves the spontaneous breaking of Poincare symmetry, it 
could persist in the double scaling limit of this model in d = 2 |193[ 1191] , 
where we also studied a NC linear a- model with 3 flavours |191] . 

Our next focus in this framework was NC QED in d = 4 |194j (with 
two NC directions). We measured the photon dispersion relation, which 
is distorted compared to the commutative world. A non-perturbative re- 
sult for this distortion is of phenomenological interest: it was suggested 
to be used to extract bounds on the non-commutativity parameter from 
the confrontation with experimental data, obtained in particular from 
cosmic gamma rays |195] . For instance, the Gamma-ray Large Area 
Space Telescope (GLAST) project (see e.g. Ref. |196j ) was scheduled to 
be launched in September 2007 and to monitor gamma rays in the range 
from 20 MeV to 1 TeV. A relative delay for photons emitted simultane- 
ously by Gamma Ray Bursts or blazar flares could hint at a NC geometry 
— or the absence of a detectable delay could impose a narrow bound on 
the possible extent of non-commutativity. 

On the other hand, perturbation theory showed that the NC photon 
is asymptotically free [197] , but it suffers from IR instability [198] , unless 
one proceeds to its supersymmetric version. If one is not willing to put 
SUSY as another hypothetical concept on top, this property seems to 
be a true disaster for the scenario of a NC world. However, we observed 
non-perturbatively at moderate coupling a phase with broken translation 
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symmetry, which appears to be IR stable |194] . Several observables fol- 
low a scaling law, which leads to this broken phase in the simultaneous 
UV and IR limit at constant non-commutativity. Thus the double scaling 
limit leads to a phase which was not visible in perturbation theory, and 
which could accommodate the NC photon. 

We applied similar methods to study the "fuzzy sphere" formulation 
|199] of the 3d theory [200] . Further details about that study can 
be found in the Ph.D. thesis by J. Medina, quoted in Ref. [37j (that 
reference also collects further simulations results on the fuzzy sphere). 
The fuzzy sphere regularises quantum field theories with a finite set of 
degrees of freedom. Hence the hope is that it might be useful for Monte 
Carlo simulations, as an alternative to the lattice. A potential virtue is 
that it preserves a number of continuum symmetries, which are broken 
by the standard lattice formulations. In the X(f)^ model this formulation 
turned out to be tedious to simulate, but it may be in business for issues 
where the lattice faces severe difficulties, like SUSY (SDIEIIES]- Therefore 
we now proceeded to the simulation of a fuzzy sphere formulation of a 
2d Wess-Zumino type model j20l| . (Related approaches were already 
presented in Refs. [Ml [65].) This discretisation is indeed numerically 
tractable, which enabled us to explore the phase diagram [66] with a 
moderate computational effort. Again the ordered regime is divided into 
a phases of a uniform and non-uniform order, as the coefficients of the 
polarisation tensors show. 
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